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Abstract 

On  the  basis  of  a  critical  survey  of  earlier   papers  it  ie  eliowri  that,  the 
lesvilts  obLaiiied  in  these  papers  are  insufficiently  justified.     A  rigorous 
solution  of  the  problem  is  given.     It  is  also  shown  that   in  the  limiting  case 
of  very  high  frequencies  the   solution  takes  on  a  furni  wliich  corresponds  to  the 
approximation  yielded  by  geometrical  optics. 

At  first  sight  the   problem  at  hand  appears  to  be  almost  trivial,  since  the 
variaDles  in  the  fundamental  equation  can  be   separated.     However,  it  is  not  at 
all  simple  to  choose  special  solutions  which  can  be  used  to  construct  a  solution 
Ox  the   pix)bltira  satisfying  all  requirements  including  the  radiation  condition  at 
infinity  and   analyticity  at  the  focus  of  the  cylinder;     nor  is  it  a  siinple  matter 
to  construct  the   solution  itself.     In  particiilar,   to   the   best  of  our  knowledge, 
the  papers    [l-U2  devoted  to  the  wave   equation  for  a  paracolic  domain  do  not  con- 
tain  a  suitably  justified  solution  of  the  problem  at  hanJ    .     In  the  present  paper 
this   gap  has  been  filled.     It  is  also  shown  that  in  the  limiting  case  of  very 
liLgh  rrequeucies  the  solution  taket.  uu  a  foim  wnich  corresponds  to  the  approximation 
yielded  by  geometrical  optics. 


1,       Fonnulation  of  the  problem 

We  cOTi&ider   the  protleu   of  ie flection  of  electromagnetic  waves  from  a  conduct- 
ing screen  wnich  has  the   form  of  a  paiaoolic  (yliuder.     We   assume  that  the  source 
of  vibrations  is    a  line   source  located  in  the  interior  of  the  cylinder  on  the   focal 
liiie  oi  the  cylinuer  and  that  tic   intensity  of  the  current  in  the  soiree  is 
I  «  I  e''''"^,  where  I,  the  aiaplitude  of  the  current,  is  constant  and  co  is  the  circular 


frequency  • 


Por  details  see  Chapter  1  of  our  paper  [5]. 
The  case  of  an  arbitrary  location  of  the  source  is  considered  in  section  8. 
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We   choose   a   syet/em  of  Caj-tesian   coordinates   (x,y,2)   so  that  the   x-axis 
lies  in  the  plane  of   symmetry   of  the   parabolic   cylinder  and  the   ^-ax;is     coin- 
cides with  thfc»   current.      In  addition   to   the  Cartesian  coordinatea  we   shall  ute 
parabolic  cooi-dinatee   (a,p)   del'iiied  by    the  relation 

(i.l)  X  ♦  iy  -  l/2(a  +  i^)^,  -oo  <  a  <  +00,     0  <  ^  <   00, 

which  implies  that 

(i.2)  X     -     l/2(a'  -  p^),  y     -     op. 

The  equation  of  the  svirface  of  the  cylinder  is  then  seen  to  be  p  -  p^  (Fig.  1), 

In  the  case  under  consideration  the  elbctric  field  vector  K  has  only   one 
component  E     of   the  loriii 

(1.3)  E       -     £(x,y)e  a     Ee  , 

Separating  out  the  field  of  the   source  we  put 

(l.U)  E     -     E^  ♦  u, 

whore 

(1.5)  E^  -  -nwlc"^  H^^\kr)i       r  =  /?7    -  l/2(a^>p^), 

(2) 
Wiere  k  =  w/c,  c  =  the  velocity  of  light  in  empty  space,  H^     (z)  is  a  Haiikel 

function  of  the  second  kind  and  u  is  the  secondar/  electric  field. 

Since  the  electrical  field  B  must  satisfy  the  Kelxnholtz  equation 

A  F  +  k     E     -^     0  as  well  as  the  radiation  condition  and  since  it  must  vanish 

on  the  surface  of  the  conuucbiag  screen,  we  can  give  th«  following  mauheniati  val 

lomiulation  of  the  problem  under  consideration:     find  a  solution  E(a,,8)  of  the 

equation 


2a  - 


Figure  1 
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(1.6)  {dh/^a)   *   (a^£/dp^)   +  k^(a^  *  p^)E     -     0. 

which  is  regular  in  the  interior  of   the  parabolic   cylinder    Qdth  the  exception 
of  the  line  o  ■  p  -  0  in  the  vicinity  of  which 

2  .2 


(1 


.7)  E     .     -^    h1^)(.14^). 


(where  u  is  a  regvilar  functioaTJ ,  and  eatisfiea  the  boundary  condition 

(1.8)  El  -0 

and  the  radiation  condition 

(1.9)  El  -0(|al-^/2),Ag*  ikE^  .0(|ar^/2) 

llaKoo  V^^  yial^oo 

(uniformly  in  p). 

In  the  sequel  (cf.   section  6)  we  shall  show  that  in  this  fomulaticn  the 
problem  has  a  unique  solution. 

2.       Construction  of  a  formal  solution  of  the  problem 

To  achieve  separatiT'n  of  variable*  in  the  wave  equation 

(2.1)  oVaa^)  ♦  (dV^p^)  *  k^(a^  ♦  P^)u     -     0 
we  put 

(2.2)  u     -     A(o)BO), 

Then  we  obtain  for  the  factors  A  and  B  the  ordinary  differential  equations 

The  relations   (1.9)  follow  in  the  usual  manner  from  the  condition  that  in  a  unit 
of  time  the   same  amount  of  energy-  radiates  to  the  outside  across  every   section 
a     ■     const,     of  the  guide  and  that  no  waves  are  reflected  back  from  infinitiy. 


-  u  - 

(2.3)  a"    +   (kV  ♦  \)A  -  0}  b"    +   (kV  -  ^)B  -  0, 

where  X  Is  the  separation  parameter. 

In  the  first  of  these  equations  we  make  use  of  the  substitution 

(2.U)  A(a)  -  exp  J-  iko^/zl  v(c)j  «  -  y^  a, 


and  we  put 

(2.5)  X    -    ik(2v*  1), 

this  yields  the  following  eqiiation  for  the  function  v(4) 

(2.6)  ▼"    -  2€v'  +  2vv    =    0. 

The  general  integral  of  equation  (2.6)  is 

(2.7)  t(o-«f('-^.|.c^\   *m(^.}.A. 

where  M  and  N  are  constants  and 

(2.8)  F(a,b,z)     .     EM    £    IVirO       z^ 

r(a)     S=-o    p(b*n)       "^ 

is  the  confluent  hypergeometric  function. 

Since  in  our  case  the  solution  must  be   am  even  function  of  the  variable  o 
we  must  have  N  '  0,  so  that 

(2.9)  A(o)  •  A^^ha)  -  exp  J-  ika^/2  I  F(-v/2,l/2,ika^) 

Here  we  may  restrict  ourselves  to  the  interval  0  <  o  <oo  which  is  what  is  done 
in  the  sequel. 


cf.,  for  instance,    [6,7]. 
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If  we  make  use  of  the  asyaptotic  represantation  for  the  function  A   (a) 
for  a  -**oo  ,  namely 

(2.10)    iS^ha)  I   ~0(v,a)exp  jika^/2|.  a"^^*^^  +  T(v,a)exp  ^-ika^/2 

a->c» 


where  0  and  T  are  functions  which  tend  to  finite  limits  ae  a  -i>co  and     if  we 
require  that  the  special  solutions  which  we  are      seeking  be  in  the  nature   of 
wa^ies  going  out  to  infinity  with  aaiplitudes  tending  to  zero,  then  we  obtain 
for  the   paraneter  v  the  condition 


(2.U)  -  1/2     <     Re  -(v    y    <  0 


For  a  chosen  value  of  \  the  general  integral  of  the  second  of  the  two 
equations  (2.3)  is 

(2.12)         B  -  B^(p)  -  exp  J-ikp2/2l  \pt(^  ,  ^  ,   ikp^)  *  Q /Ik  pF  (^1  *  J,   ^  ,ikp^)], 


where  P  and  Q  are  constants. 

In  view  of  the  bcundedress  of  grad  u  at  the  focus  it  follows  that  Q  »  0  ,   so 
that 

(2.13)  B^(?)  .  b[^\^)  '  exp  J-ikp^/2l    f(^,\,   ikpM    . 

Thus  we  arri've   at  the  special  solutions  of  equation  (2.1) 
(2.II4)  u^  .  A^^\a)B|,^^(p)j       -  1/2  <  R  J   V   I    <  0, 


and  a  possible  soluticn  of  our  problem  is  given  by  the  complex  integral 


See  Supplement  1,  formula   (D,l)  as  veil  as    [8] 

«*  This  result  is  obtained  by  substituting  the  special  solutions  under  consideration 
in  the  equations  (3.5)   in  the   sequel. 
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-6  +ioo 
(2.15)  u(a,p)     -      j  C(v)A^^^a)B^^^O)dVj         0  <  6  <  1/2. 

■Lb 'ioo 

Ihs  unknown  function  C(v)  must  be  deteiinined  from  the  boundary  condition  (1.8), 

i.e.,  from  the  equality 

-6 +100  ^22 

need  ,(2)  /     °     *  P. 


.16)      I  C(v)A^^\a)B^^^O^)dv 


(2 

^6 -ioo 


C3J 


Thus  the  problem  is  that  of  representing  the  field  of  the  source  as  an  integral 

with  resi»ct  to  the  functions  a},     (a).     Such  a  representation  can  be  obtained  by 

making  use  of  Magnus*  deconpositlon  theorem 

-5+ioo 
(2.17)  f(a)     -     I  F(v)A^^^(a)dv; 

-6 -ioo 


pL     p(.  v^  rC^),^--/^  j    f  (a)A^^^a)da       (0  <  6  <  1), 


2n' 


from  which  it  follows  that 


Making  use  of  the  fact  that 


^^^"  '^'  .1/2  i«/Un(3) 


(2.19)      (  H(2)(^,5!j£Wl)(a)da-(nk 
where 


)--/'  e-"'-*  b;^'(p), 


where 


B^^h^)  -  exp|Iik^2  /sj^  p  F(l  ^  ^  ,  ^  ,  ikp2) 
is  an  odd  solujiion  of  the  second  of  the  equations   (2,3), 
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(2.20)  bI^H^)     '      (iK^^)--^/'^    r(^)  W     ^       ^       -^(ikp^) 

(W,      (z)  is  Whittaker'E  function'-  -^ )  we  ootain  the  final   form  of  the  solution  of 
our  problem 

(,.21)  u(a,^)   .   ^1-7  r(-  ^)r(^)e^^'^/^  !^B^^)(,)Ai^^a)d. 

(0  <  5  <  1/2). 

It  is  ciear  that   this  formal  solution  requires  detailed  verification.     This  veri- 
fication is  found  in  section  5. 


3,       Investigation  of  the  integral  representing  the  secoitdary  field  and  proof  of 
the  regularity  of  the  obtained  solution 

1,     The   functions  which  appeai   in  the  integrand  in  (2.21)  are  entire  functions 
and  as   such  are  continuous  functions  of  tljti  variable  v,     Gince  the  function 
B       (^    )  does  not  vanish  along  the  path  of  integration   (see  supplement  2),  it 
follows  that  the  integral  exists  for  every  choice  of  finite  limits  of  integration. 


shall  make  usa  oi  th«  asymptot-ic  repr-es«iitations  for  the  functions  k       (a),  B       (p), 
and  B       (p)  for  fixed  a  and  p   (A  <  A,  0  <  |^  <   oo),    iv|  ->oo,  R^v^'  >  -6  > 


To  prove   convergence  of  the  integral  for  infinite  limits  of  integraticn  we 

nctic 

,.f 

Putting  V  =  -5  +  if,  we  have  for  T  -^--oo: 

(3.1)        |Ai^^a)|  =0(l)i    \Bl^h^)\   -0(«^^  >^^TT4  ]i 


|Bp^?)| 


Since 


(Ylrf^/^  expl-vnTTHPn 


cf.  formulas  (D,3),   (D,U)   and  (D,$). 


I  nc-^JTc^).'""^  I  -0(1^ 


-v.^ 


it  follows  thai 


-6-iT  ,(3), 


which  implies  the  uniiorxa  convergence  of  the  integral  (2.21)  with  respect  to  the 
lower  integral. 

In  Just  the   same  way  we  find  lor    'C  -^  +00: 

(3.3)  |A^^^a)|   -  0|e^^    /^Ffn    }  ]i 

\B^^\^)\  -0(f  ^/^); 


whence 


(3.U) 


,b(3)(p),  .o(4i)i  In-  J)r(^)e'"'/'  -ar-^^'e'^' ) 

<.6+iT  ^v     ^^0' 

,00 

,.    I       -nr     /2kr       df 
<  M   I     e  e*^  -^ 


so  that  the  integral  (2.21)  converges  uniformly  also  with  respect  to  the  upper 
limit . 

The  uniform  convergence  of  the  integral  (2.21)  implies  that  u  is  a  continuous 
function  of  the  variables  a   and  ji  in  the  domain  under  consideration. 


2.     To  pi  ovfc   anal>ticity  or  tho  lunction  u  we  must  still  prcrve  the   existence 
and  continuity  of   the  first  two  derivatives  of  u  with  respect  to  x  and  y.     For 
this  it   6U If ices  to  prove  the  uniform  convergence  of   the  integrals  obtained  after 
differentiation  under  the  integral  sign. 

The  estimates 

lA^^^'Ca)!   .OCI-CDi    |Bi^^'(|i)|   -0(irl«^P  jv^^n^l^l     )   -^^     t— 00 

imply  the  uniform  convergence  with  respect  to  tl:ie  lower  limit  of  the  integrals  ob- 
tairfid  from  (2.21)  by  dif f exentiotiiig,  under   Lhe  ini/fcgial   faign  witli  respect  to  a  and 
p.     In  exactly  the   sarae  vay,   making  use  of  the  estimatee 

I         1 

|aJ,^''(o)|    -OCt-eip  jy^FTa}-);    |BJ,^*'(|i)|  .0(t)       as     T -^  ♦CD, 

we  can  prove  the  uniform  convergence  of  the  integxalt  in  question  witli  respect  to 
the  upper  limit.  Analogous  reasoning  proves  the  existence  and  continuity  of  de- 
rivatives of  any  order  of  u  with  respect  to  a  and  p. 

We  now  consider  derivatives  of  u  witii  respect  to  x  and  y.     These  derivatives 
are   given  by   the   relations 

,        .  du  g         du  p       8u         8u  a       3u  p       du 

a  +p  a  ^p       ^  a  +^       ^       a  ■'p 

liicse  foriiiulas  sliow  tha^  3u/3x  and  du/dy  are  continucws   in  the   domain  under  consider- 
ation with  the  possible  exception  of  Uie  point  a  =  |3  =  0 .     To  investigate   the  be- 
havior cf  du/ax  and  Su/dy  at  a  -  jB  -  0  we   put 

(3,6)  a  =  r  cot  y;       p  =  ^   sii'  ^         (0  <  w  <  i.) 

and  consider  the  expression 
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(3.8)  7^^(r)  =  i^os  9  B^^\r&iii&)A^^^'(rcose)-£in  &  A[^\rc08e)B^-'-^ '(rsine)! 

Since   for  r  -i>  0     /'A^)   tends  tc   a  finite  limit,  nanely,   to  -ik(2v+l),   and 
since  the  estimate 

valid  in  the  domain  r  <_  r  .iicpliws  the  uniform  convergence  or  the  integral  (3.7), 
we  may  go  over  to  the   limit  under  the  integral  sign.     Thus 

lim     ^  .  -  Jhi^  r(.  -)  r(i^)ei"^/2(2v.l)  ZJl^  dv 


Since  tlds  integi al  is  independent  ol  ©,  we  may  conclude  that  8u/dx  is  continuous 

in  the  neighborhood  of  a  -  3  -  0. 

Similar  computations  show  that  lim      ^  =  0  wiiich  proves  the  continuity  of 

r-oo     "^ 
Bu/8y  in  the  nei^borhood  of  the  focus.     Using  similar  metnods  we  can  prove   the 

continuity  of  the  second  derivatives  of  u  with  respect  to  x  and  y. 


U,       Varioas  representations  of  the  solution 

Before  we  can  show  that  u(a,0)  is  actually  the   solution  to  our  problem  we  must 
represent  it   in  a  number  cf  ways. 

1,     To  obtain  the  first  of  these  repres  ntations  we   complete   the  path  of  inte- 
gration in  formula  (2,21)  by  adding  to  it  a  circular  aic  of  ladius  R  -p-oo  and 
center  v  =  -  1/2  and  utilize  the   residue  theorem.     Asymptotic  expressions  for 
functions  of  a  parabolic  cylinder     show  that  with 

*  CI.  (d,3)  -  (D,5K 
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(U.i)  a     <     2^^  -  p 

the  value  of  the  integral  over  the  circular  arc  tends  invariably  to  zero. 

If  the   inequality    (U.I)   holda,   then,   computing  the  residues  at  the  poles 
V  =  2n,  we  find  the   following  repreeentation  for  ui 

A  sindliir    Liane formation  of  the   contour  integral  (2.16)-(2 .15)  for  the  primary 
field  yields  the  repi esentition 

2\         21       ^       fin*  I) 
Pin  *1) 


valid  up.der  assiunption 

(Ij.h)  a'<     p      . 

We  observe  that  there  is  a  simple  c<3nnection  between  the  functions  A^     (a) 

and  Herxnite  polynomials,  namely 

ih.'^)  A^^a)    =   (-1)"  ^     exp  |-ika2/2]     H^   (  >^k  a) 

inus  at   least  for  a  <  p  our  solution  can  be  represented  as  a  power  series  in 
Hermite  polynomials  H^   (  y5T  a).     As  for  the  domain  a  <  p,   direct  investigation  of 
the  convergence  of  the  sezies   (U.3)  shows  that  for  a  >  p  this   series  diverges. 
Similarly  the  series  (ii.2)  diverges     for  a  >  2pp-p.     We  may  thus  claim  that  the 
total  field  can  be   repretented  iu  a  series  in  Heri.iite  polyacniialfc  only  in  the 
domain  defined  by  the  inequality  a  <  |3.     In  other  words,  there  exists  no  single 
analytical  representation  of  the  required  solution  in  the  form  of  a  series  in 
Hermite   polynomials     • 

2,      ]jr\  the  paper   []J  Magnus  gave  the  following  solution  of  the  problem  under 
consideration: 


Convergence  of  these  series  was  investigated  using  asymptotic  expressions  for 
Hermite  polynondals  for  n  -i>oo, 
*<  It  appears  that  this  fact  has  not  been  properly  noted  in  the   literature. 
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.-e*i-    3(3)(,  J 


-6-ioo       \     ^P- 


(U.6)       u  -  -^  \.     °     Bi^h^)ki^\  |a()dv       (0  <  5  <  1/2) 


2-^  4-ioo    b;-^(p,)  ^       ^ 


where 


(U.7)  A^^ha)     '     B^J^^O). 

Since  Magnus  does  not  prove   Uiat   the   solution  given  by  him  is  regular  everywhere 
in  the  interior  of  the  domain  under  consideration  it  is  impossible,  without 
further  investigation,  to  conclude  that  his  solution  is  correct.     The  need  for 
such  investigation  is  enliaxiced  by  the  fact  that  the  special  solutions 
A^     (ia|)B^     (p)  chosen  by  Magnus  are  not  regular  in  the  neighborhood   of  the 
focus  a  =  ^  =0. 

That  the  solution   (U.6)  given  by  Magnus  coincides  with  our  solution  (2.21)  in 
the  domain  in  which  both  solutions  are  regular  can  be  showi  by  adding  a  circular  arc 
of  radius  K  ->oo  on  the  right  side  of  the  path  of  integration  in  (U.6)   (which  is 
permissible  in  view  of  the  fact  that  a  <  2p  -[j)  and  by  computing  the  residues  at 
the  pole 6  V  -  2n»     These  computations  lead  to  the  formula  (U»2)»     In  view  of  the 
well  known  properties  of  the  solutions  of  elliptic  equations  we  msy  thus  conclude 
that  the  solutions  (U.6)  and  (U.2)  coincide  in  their  common  domain  of  regularity. 
Finally  ve   observe  that  the  field  of  the  source  can  also  be  represented  ab  an 

integral  in  the  functions  A^^^(|a|): 

,-6+i  00 
(U.3)  ii\'-'lk^-^\'^\  Bp^p)Ap^(|a|)dv       (0<6<1), 


"i"e-^)-*L 


and  this  representation  is  valid  in  the  whole  domain  with  the  exception  of  the 
line*  8-0 


cf.    [U] 
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3«     On  the  basis  of  t\ie  results  in  section  2  the   complete  solution  E  can  be 
ropr^faeiited  by  weans  o£  tlia  integral 


r^*ioo  .(3)/,    |N 

\~J^)f.    ..(?).. ^_o(2),.    ..(l)/..Tv  Jlllildv     (p-0), 
-6-ioo  "  ""v     ^f^o^ 


-0*100  U- 


which  permits  us  t:>  obtain  another  form  ol  the   solution  of  great  use   in  the  sequel. 

We  now  complete  the  path  of  integration  by  addition  on  tht!   left  of   s  ciicular 
arc  of  radius  R  -.>oo  and  we  use  the  residue  theorem.     Since  the   singularities  of 
the  integrand  are   the  roots  of  the  equdtion 

(U.IO)  si^^(^o^     "     °  (n  =  1,  2,   3,    ...), 

n 

located   (see  supplement  2)  on  the  line  Re  |  v  >      •  -1/2,   and  since  Im    ',  "^i  \  <  1">  |  v^  \ 
<   ...,  anJ,   finally,   since  the  iiitegi  al  along  the  arc  tends  to  zero,   it  lollovts  that 


(u.u)  E-i^^  f    T-TD ~ ^?^(l<^l) 


r—  n  n 


Computation  of  the  derivative  which  appears  iii  the  denominator  leads  to   uie 

followxng  representation  of  E  in  a  series  in  the  functions B^     (p)» 

n 

(U.12)  E(a,ft)   -  -  ^-^^    r    B^^^(p)A^^^(|a|)/    (  %^^^^(^)dp    . 

n^       an  In 

This  form  of  the  solution  will  be  utilikied  in  the  next  section  \rfien  we  prove  that 
the  solution  of  oui    problem  satisfies  the  boundary   condition  and  the  radiation 
principle. 

We  wish  to  note  that  the  representation  (U.12)  can  be  obtained  directly  by 
iiilegictiiig  the  Wpve   equation 

Here  we  made  use  of  the  relation: 


Hi  - 


(U.13)  {b\/ba)  +   (8^E/dp^)  +  k^(a^+p^)£  -  Unicoc*^(a^+p^)j   , 


\inder  the  absumptioft  t2iat  thd  intensity  of  the  current  is  vmiforiQlir   distributed 
over  a  small  neighborhood  of  the  focus: 

1/2  {a^^<^^)a^^^         for    |a|   <  a^^,     O  <  p  <  ^^ 
j(a,p)  -^ 

0  for  the  remaining  values  of  a  and  ^ 

with  o^  and  p     tending  to  zero. 

5,       Test  of  the  formal  solution 

The  results  of  section  3  imply  that  the  function  u  constructed  above  and 
defined  by  formula  (2,21)  is  a  regular  solutiun  of  Heliaholtz's  equation  at  every 
point  of  the  domain  under  consideration.     As  regards  the  complete  solution  E  it 
is  tlie  sum  of  the  field  of  a  line  source  througli  tlje  focus  and  the  function  u  and 
as  such  it  obviously  satisfies  the  Helmholtz  equation  at  all  points  of  the  domain 
other  than  the   focus  where  it  has  a  singularity  of  the  required  type. 

It  therefore  remains  to  check  if  the  soluticn  E(a,p)  satisfies  the  boundary 
condition  and  the  radiation  condition.     To  do  this  it  is  best  to  use  the  form 
(U.12)  of  the  solution. 


J^^\p)  satisfy  the  condition  B^,-"-^ 
a  n 

ly  in  order  to  check  if  the  boundary  condition  E(a,^   )  =  u  is  satisfied  it  suffices 


i.     Tne  eigenfTinctions  B       (p)  satisfy  the  condition  B^     \^q)  "  ^*     Consequent- 
n  n 


o' 


Using  the   asymptotic  expressions  for  the  functions  B^  ' {^)  and  A^-^^a)   for 

V     -  -  1/2  +  i  t       (  f      ->oo): 
n  n  n 

(5.1)  B^^^^^)  c^  cos  /5i:rr;;  PJ     k^^ha)  ^    2  /2ni/  T^  exp  J-  V^k^  al  , 
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and  the  formula  (D„2) 


(5.2) 


1         ,2n-l     ^2 
7     (-5—  «)    » 


we  find  that  the  terme  of  the  series  (ii.l2)  are  of  the  order  of 


n-e.  -. 


whicli  implies  the  unifora  convergence  of  the  ssrlee  in  the  indicated  domain. 
This  shows  that  our  solution  satisfies  the  prescribed  boundary  caidition, 

2.     To  see  that  our  solution  satisfies  the  radiation  condition  (1.9)  it 
suffices  to   show  that  the  following  two  ecjialities  hold 


(5.3)  Ul       __  -r)(»"^^^):        I  i|i*iku 


l"l..>oo  -0(°-^^')i 


a  da 


-  0(a-l'''). 


From  the  estiwatfes'  (n  -t-oo) 


(5.U) 


B<»0) 


'  n 


-OU).        |A"'(<.)|-<.-^^^0(n-5/^). 


we  conclude  that  the  first  part  of  the  radiation  condition  is  satisfied. 

W.fferentiating  the  series  (U.12)  and  using  the  estimate   (deriwd  in  supple- 
■ant  3) 


(5.5)  |i»i"'(«)^1^4''(«)|„_--5^^0(n-5^^) 

'  n  n  ' 

vTliich  is  valid  for  arbitrary  a  we  find  that 


The  second  estimate  (5.1),  derived  in  supplement  3,  is  uniform  vn.th  respect  to  a 
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'a-i>oo 


which  proves  that  the  second  part  of  the  radiation  condition  is  satisfied. 

Thus  the  formal  solution  of  our  problem  obtained  in  section  2  satisfies  all 
conditions  specified  in  section  1  where  the  problem  was  mathematLcaUy  formulated. 


6»       The  uniqueness  theorem 

A  solutiOii  of  the  Heluholti.  equation 

(6.1)  /^  V  +  k^v     »    0, 

which  is  regular  in  the  domain  ( D)  bounded  by  a  branch  of  the  parabola  P  ■  p., 
and  the  axis  of  abscissas  and  satisfies  the  boundary  conditions 


(6.2)      v|p.p^  -  0,      dv/8y|y_^  -  0 

•«■ 
and  the  radiation  conditions 

vanishes  identically. 

For  p:'oof  we  consider  the  domain  (C)  bounded  by  the  x-axis,  by  a  branch 
of  the  parabola  p  =  p     and  by  a  branch  of  the  paraoola  a  =  A  ortnogonal  to  the 
parabola  p  ■  p   ♦     Multiplying  equation   (6.1)  by  the  conjugate  function  v  and 
the  conjugate  equation  by  v  and  integrating  the  difference  of  the  two  ecjiations 
over   (D»)  we  find,   using  Green's  theorem  and  the  conditions  (6.2),  that 


(6.U)  J 


Using  a  consequence  of  the  radiation  condition 


Uniformly  vdth  respect  to  the  variable  p< 


-  17  - 


(6.5)  I   U|I-  V  |I  j  .  -  2ikivi^  .  vo(a-'^/^)  -  vo(a-^/'^) 

and  noting  that,  I'or  o  -*-oo  the  products  a       v  and  a       v  remain  bounded  we  obtain 
from  (6.U)   tlie  lollowing  equality: 

(6.6)  liB  a|v|^  dp     =     U. 
o->oo      J^ 

To  show  that   (6.6)  implies  that  v  •  0  we  consider  the  following  special  solutions 
of  t.iie   rieimnolx,z  equation: 

n  n 

where  the  v     are  roots  of   the  equation  B^     (Po^  ^  ^*     Ti'**^«  solutions  aae  regular 


X-axis     and  satisfy  the  baundar>   coi^Jitions 

and  the  radiation  condition   (6.3). 

Applying  Green's  formula  to  the  pair  of  functions  v,  u     we  obtaiii  the  relation 


J     ;    n  da 


(6.y)  ('Yu.  g  -  V  ^A  dp  .    rYu„  i  -  '  ^  1  ^.,     "P' 


This  relation  implies  that  the  value  of  the  integral 


rf„E 


du 

"o 

is  independent  of  a.     In  view  of  the   fact  that  this   integral  tends  to  zero  as 
a  -*•  00   (   a  consequence  of  the   radiation  condition)  we  conclude  tiiat 
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(6.10)  I        (   u     g-v^    Wp     ^     0         xor  alia. 


(-(:■ 


SutJsiituting  for  u     its  value  aiiu  putting 


(6.n)  v^     -      ]     V  B^^>dp  /[       B^,^^'  dp  (n  -  1,   2,   3,    ...). 


n 


whenee 


n 

Niiere  c     is  a  oonsbaxit, 
n 

(l)  * 

Giace  tfte   grstem  of  fvinctions     B       (p)  is  closed     ,  i.e., 

n 


(6.13)  Z.    K\"    I     ^i^^  ^P     •     I       I'l^  ^P' 


2 

we  conclude  fi-cm  (,6.o)   Uial     lia       al'^'nl     "^  *-*•     ^  view  of   (6.12)   tiiis  meaiiB  that 
a-i'OO 


(6.1U)  Ic    i^     lim       a\f^l/'\^     -    0. 


(3)  ♦* 

Now  the  asymptotic  representation  of  the  functions  A*"   (a)  fur  a  -^oo   show£  that 

(3)    2  ^ 

the  piouuct  o|A^,  (a)|  does  not  tend  to  zero  for  a  -»>oo.  Hence  c  =  0,  i.e., 

n 
(6.15)      \  ■  0;     n  «  1,  2,  3,  ...  . 

In  view  of  the  closedness  condition   (6.13)     v  -  0  for  all  a  >  0  and  0  <  p  <  p   , 

This  follows  from  the  general  theory  of  eigenf unctions  associated  with  second 
order  differential  equations 
•Mt  cr.   (D^2) 


-  ly  - 

which  is  what  wo  wished  to  prove. 

7 .       InvBstigation  of  the   aolutlon  for  large  valxas  of  zhv  >.ove   numbtii    k,     Itanfei- 

tlon  to   geometrical  optics. 

We  now  ia^stigdte   the  behavior  of  the  solution  ootaii^d  in  beclxon  2   aS 
k  -->oo  widch  should  correspond  to  the  case  of  geometrical  optics.     To  study  the 
behavior  of  t-he   secunuaxy   field  as  k  -i>oo  we   make  use  of   the  following   aeymptjotic 
formvilas 


'i,"<s)U  ^co  ~n^)(V)-<'"*^V''"' ^'. 


ry 


whicii  are  obtained  from  (D,l)   and  (ll,2)  by  replacing  v  b>  -v-1  and  tsy  noting 
that  A^(a)   -  B_^_^ia) , 

Substituting  (7.1)  into   (2.21)  yields  the  relation 

«^<^'^Hk  ^00     ~     -f-T exp^^(p-.p^A     i^   i 


-6+100  ^ 


To  compute  the  integral  I^   we  make  use  of  the  integral  repiesentation   (iLoy,  widcL 
when  put  ill  (7.2)  yiflds: 


-  20 


We  now  make  the  substitution  v  =  2n^l  in  Ihe  complex  integral.     If  ve  put 
26'   "  1-5,  we  find  that 

2p2  ^'+ioo  k 

h  ■  li  «''(-  r  - «'  JTT  J      r(.)(-ike-»  -I )-"  d., 

'^  6'-iQ0  P 

whence 


M 


(7.U)  I.  -  2  yr    (^)  exp  J-  ^  -  9  .  ik  ^  e-2»l 


Substitutiiig   (7.U)  into   (7.3)   and  effecting  the  change  of  vai-iabl©  e"**  *  x  in 
the  Integral  I.  we  find  that 


(-)  "iw..  ^-  ^(.^-^.^  |eJ- ^-^^ -. ;44^ 


As  k  -.>oo  the  principal  contribution  in  tiie  latter  integral  comes  from  the  nei^- 
borhood  of  the  point  x  =  0.     Hence 


(-^.6)  u|^^^  -  i-  2l(nk)^/2g  ,(,^/^^)2J"^     expJ^  -  ikp^.  ^  (^V) 


If  we  now  revert  to  Cartesian  coordinates  and  note  the  relation  (of.  i-'ig.  2) 


5 » +1  CO 


y5»-i: 


--^  -     "■       '  PM  x-^   dp 


-  2ea 


Fjcrure  2 
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we  can  rewrite  formula  (7.6)  as  follows: 


(7.7)       u|,_  c  ^   ^-^-^ 


From  Llwj  laiter  relation  it  is  apparent  that  the  exponent  s  widen  appeals  in  tiie 
phase  idctoi'  reprefcents   Lhe  path  0AM  wiiich  would  be  tiaversed  by   a  ra>  emanating 
from  the  origin  U  and  pabsiag  through  M  as  a  resvilt  of  reflection  from  the  surface 
of  the  parabolic  cylinder,   such  reflection  obeying  the  laws  of  geonetiical  optics, 
■^t  is  possible  to   show  that  the  amplitude  of  the  factor  in  formula  (7.7)  also 
agrees  with  the  appioxLmation  yielded  by  geometrical  optics.     Indeed,  if  we   compute 
the  mean  density  p  over  a  period  of   lhe  aieigy  flow  at  some  point  M  ; 

(7.5)  p     -     (c/l6n)(u*H  -  uhJ), 

and  if  we  make  ute  of  the  asymptotic  representation  for  H     (wriich  is  identical 
with   (7.7)),  we   find  Lhal 


(7.y: 


PMco^^^l'/"P^/[^*   (^/^o4 


un  the  other  hand,  we  are   led  to  tids  very  exprefcaion  by  assuirdng  the  laws  of 
geomeoiical  optics.     To  see  this  we  compute  the  electromagnetic  eneigy 

IjL  d»  -    ■    i'*^    ,     ..      ju^l    -  pidy  i  wiiich  pasefcs  in  a  unit  of  time  through  a  small 
Zn  2n|ay/d©| 

area  of  height  dy  lu  t.ne  vicinity  of  the  point  M  (cf .  Fig.  2).     Here  W  -  /I  kl  /2c 
denotes  the  average,  with  rebpect  to  time,   amount  of  encigy  radiated  Dy  a  unit 
IciigU.  oi   the  source.     Using  the  geometrical  relations 


neve  we  discard  tije   component  H     of  the  secondary  magnetic  field  since  the  latter 
is  of  the  order  of  1/k  in  coiuparisoii  with  H.  . 


tg»  •  y/x,     y^   -  2p^(x+p^/2)  we   find   that 


o2     ^    w      dy 


^  E  •  t>^]  • 


;«lu.ch  leads   Lo  the  expression   (7,9)  for  p. 

Thus  in  the  limiting  case  k  -»oo  the  solution  of  our  pioblera  poee  over  into 
the  approximation  yielded  by  geometrical  optics. 

8,       Generalization  of  the  solution  to  the  caee  of  arbitrary  location  of  the   source 

on  the  axis  of  the  cylinder. 

In  the  cabe  when  the  source  of  vibrations  is  not  located  at  the  focus  the 
method  of  solution  presented  in  section  2  leads  to  complex  computations  connected 
with  the  repi«8entation  of  the  field  of  the  source  as  an  integral  in  the  functions 
of  a  paiaDolic  cyliader.      In  cuch  eases  it  is  expedient  to  look  at  once  for  the 
complete  field. 

We  first  consider  the  case  when  trt*   source  is  located  at  an  arbitrary  point 
a=o,^=Oof  the  positive  part  of  the   axis  of  abscissas.     Replacing  the  lir^aar 
curieal  I  by  a  current  equal  to  it  but  distributed  uniformly  over  a  small  area  we 
csn  write  the   following  non -homogeneous  equation. 

(5.1)  i^.i^.  Aa'  .  ^^)E  =  Itl^  (a'  *  ?')i(a,^) 

da         Sp  c 


(0<a<oo,      0<p<  p^)| 


1/     U(a  *^   )eY  ior         o*  -e  <  a  <  a*  +e,     0  <  p  <  y 


0  otherwise 
and  the  boundary  conditionj; 

(8.2)  E|p^^=    0,      aE/a.|^^^.    C-   aE/9p|p^=    o 

ir  we  seek  the  function  E(a,e)    in  the  form 
■5+ioo 

(8.3)  E     =     Ji 


^ E^(p)r^^)n(^)e^""/'Ail)(a)dv,      0  <  6  <  1/2   , 

2n  I_5-loo 


where 


1  -<^^ 


23  - 


(8./*)  Ejp)     =     J     EAl-^^(a)da, 

0 

then,  multiplying  (8,1)  by  A^  '(a)  and  integrating*  witli  respect  to  a  fron  0  to 
OD  w©  do  tain  for  E   (p)  an  ordin.iry  differfintial  equation.     Let  E  (8   )"0,  E'(0)  =  0. 
Under  the  assumption  that  e,  t  -*"  0»   these  conditions  imply  that 

(8.5)  Z^O)  -  5fE  Bi^'(p)B<3)(pJ   -  b(3)(„b(^)(,^,  .^    . 

Substituting  (6.5)  in  (8.3)  we  obtain  the  aolution  of  our  piooieiu  • 
-6+ioo  _ 

2Ti/i;    C         J  .  l- 


-6+ioo 

I    f 

-6-xoo 
(6.6)  ,(1) 


Bv     (P)B^     Oo-U     (1),     ■■  K     ^°^"^* 


-V     (^o) 


An  analogous  procedure  can  be  used  in  the  case  when  the  source  is  locaued 
to  the  left  of  the  focus  at  an  axbitrary  point  of  the  negative  part  of  the  axis 
of  abscissas. 

In  the  preceeding  discussion  we  have  restricted  ourselves  to  the  case  when 
the  source  is  located  on  the  axis  of  the  cylinder.  A  similar  dismission  could  be 
applied  in  the  case  when  the  source  is  located  in  an  arbitrary  position  inside  the 
cylinder.  In  this  more  general  case  one  would  have  to  introduce  in  ^.H-iition  to  the 


even  f\mction6  A^  (a)  also  odd  functions  for  which  it  would  be  necessary  to  derive 
the  corresponding  representations. 

*  Here  we  make  use  of  the  radiation  condition 


2U 


Supplement 8 


1,       Some  propwrties  of  the   functions  of  the  parabolic  cylinder, 

AsyBiptotic  behavior  for  large  values  of   the   arguwci.t: 

2 


Ai^>(«) 


i„|.„  =  «""^^    m4i(,.„^)-(-)/^[..0(;;i.)} 


I  ^^  ^1  r[(v.i)/2]  L  ^  ^  J 

Asymptotic  representatione  as    |v|-^cx): 

(D^3)  A^^\a)  '  cos  /ik(2v*l)     a   [I  *  Q  (—-)]   » 


(D^U) 


[,^'(^)  -  ch/ik(2v*i)  p    [1.0/^-:^'\] 


(D,5)*  B^^^p)  -     ?     iiL  exp  ^   *  /ik(2v.l)  p 


Integral  representatione 


W  _  i^  th  0   I    -^ 


-00 

(-  1  <  Re  V  <  0)j 


Here  the  upper  sign  refers  to  the  domain  Re  I  vl  <  -  j    and  the  lower  sign  to  the 
domain  Re  j  v  |    ^  -  j-  • 

■JHt-  PlOl 

This  formula  is  obtained  from  the  well-known  integral  representation^    -• 

1 

F(a,b,z)   -   DLh) (     t*-^(l-t)^'^'V^dt}   (Reb  >  Hea  >  0) 

n(a)r(b-a)     Jq 

20 
by  the  ch&neie  of  variable   (l-t)/t  =  e 


25  - 


(V>   <^><»>-s^^-^'"^^-{r-¥lfJi-]-- 


(Rev  <  0), 


2,       Distribution  or  the   i^ioe  of  the  funclion  b^  '(p) 

lo  investigate  the  locations  of  the  roots  of  the  equation 

(D^i)  B^^h^)    =    u     (v  =  <r.  ir). 


we  mate   use  of  the  relation 

a 

Bi^^(a)Bi^^'(p)  -  ^^^\?)B^y{^)    =  2ik(l*2<r)    f   IB^^^I^  dp. 


U2<r)  I  |B^^' 


which  can  be  gotten  from  the  differential  equations  and  the  iidtial  cuiditiuns 
for  the  functions  under  consideration.     If  we  assume  tnat  for  a  given  p  the  value 
of  V  is  such   iii&t  B^^^h»)   '  0,  then  b!^\<^)   =  0.     Consequently,  Re[v^=  o"  -  -  1/2, 
i.e.,   all  roots  v     of  the  function  B^  ^(p)  are  locaUd  Oii  the  line  Tueivj^  -  1/2. 
A  more  detailed  investigation  shows  that  the  imaginary  parts  of  these  roots 
If    i   -i>co,  are  given  asymptotically  by 


<v)     ^n-,-^c^«y 


It  can  also  be  shown  that  the   numbers      X     are  bounded  from  below 


n 

3, 


3,       uerivalion  of  a  certain  estimate  for  the  function  A^(a) 

It  is  the  purpose  of  this  supplement  to  establish  a  certaxi*  estiiaate  for  the 
function  of  tlio   third  kind  A^-^^(a)  valid  for  -T  <  Im  j^vj  <  oo   (T  is  a  large  positiMO 
number)  uniformly  with  respect  to  a. 

Integrating  (D,7)  by  parts  we  find  that 


A0)(,).!!:i^(Ka^)V/2,^)i-.iV 


ni-v/2) 


"H T 


,00  .        2 


^    Jq  Ua  k  ^ka 

Estimating   the   absolute   value  of  the  latter   integxal   and   tdc  ing  iuLo   coij&idSi'ation 
the  behavior  of  the    P-function  for  large    T  we  fixid  that 

(D3I)  |A^^^a)|     «      a"^9f(r)  (6-0), 

and 

(0(i)  ^°^      -T  -  r  <  T, 

(D,2)  ^C^)     -     <^  ,    A/p 

^  pr"^-^/^)  iur  r<     r  <   ex,. 

further,  using  an  integral  representation  for  the  linear  combination 
a~  A^   (a)  +  ikA   (a)  (such  a  representation  can  be  easily  obtained  from  (D, ?)) 
we  find, proceeding  as  above ,  that 

(V>      '  I    ^'^^'^^  '   ikA^7\'a)  1   ^  a"'-^  I(T), 


where 


0(1)         lor   -T  <  ~  <  T, 
(D.U)      Y(t-)  -  ^  ^    .  ./2 

5  I  0<'^"     >  ^°^         T  <  T  <  00, 
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This  paper  is  devoted  to  the  problem  of  finding  the  electrcsnagnetic 
field  of  a  dipole  located  inside  a  parabolic  reflector  when  the  radiating 
dipole  is  perpendicular  to  the  axis  of  the  paraboloid.  This  problem  was 
considered  earlier  by  Pinney"-  -' ,  However,  from  the  mathematical  standpoint, 
Pinney's  representation  of  the  solution  in  the  form  of  series  of  Laguerre 
polynomials  is  not  sufficiently  well  justified  (of.  Section  J,   below). 

The  final  solution  of  the  problem, applicable  everywhere  in  the  interior 
of  the  paraboloid,  is  presented  in  this  paper  by  means  of  formulas  wliich  express 
the  required  field  components  in  the  form  of  complex  integrals.  The  exact 
solution  is  used,  among  others,  to  find  the  field  in  the  limiting  case  of 
geometrical  optics  when  the  wavelength  is  much  shorter  than  the  focal  distance 
of  the  paraboloid. 

1,  Formulation  of  the  problem  and  representation  of  the  solution  by  means  of 
complex  integrals 

We  wish  to  find  the  electromagnetic  field  due  to  a  dipole  of  moment  p 
located  at  the  focus  of  an  ideally  conducting  infinite  paraboloid  of  rotation 
and  perpendicular  to  its  axis.  Vfe  select  a  Cartesian  coordinate  system  whose 
z-axLs  coincides  with  the  axis  of  the  paraboloid,  whose  x-axis  coincides  with 
the  direction  of  the  dipole,  and  whose  origin  coincides  with  the  focus  0  (Fig.  l), 


Figure  1 
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It  follows  from  Maxwell's  equations  that  each  of  the  six  components  of 
the  electromagnetic   field  is  a  solution  of  Helmholtz's  equation.      If  we  put 

E     -  u  ♦  w  cos  2  0, 

(3.1)  E     o  w  sin  2  0, 

E     =  V  cos  ^ , 

where  u  -  u(a,ft),  v  =  v(a,p),  w  -  w(a,S),  then,  using  paraboloid  coordinates 
(a,p,^),  we  are  led  to  the  following  expressions  for  the  unknown  functions' ' 
u,  V,  and  w  . 

where  k  -  -  =  the  wave  mamber. 
c 

In  view  of  the  fact  that  on  the  surface  of  an  ideal  conductor  the 
tangential  components  of  the  electromagnetic  field  vanish,  we  have,  in  the  case 

The  nature  of  the  dependence  of  the  required  quantities  on  the  angle  0  can  be 
predicted  from  the  expressions  for  the  components  of  the  field  of  the  dipole 

(cf.  Ci,5]). 

■K-K- 

Note  that  in  a  paraboloid  system  of  coordinates  (a,p,0)  defined  by  the  rela- 
tions: X  =  a  p  cos  0,  y  -  o  p  sin  0,  z  -  i  (a  -  p  ),  where  0  <  a  <oo,  0  <  B  <c 
the  Helmholtz  equation  for  the  i-th  component  of  the  field  takes  the  form 


a  So   "35-   p  ep  '  ep 


P^/   302 


under  consideration,   the  following  two  boundary  conditions: 

irtiere  p  =  p     is  the  equation  of  the  paraboloid  and  the  interior  of  the  paraboloid 
is  described  by  0  <  a  <  oo,  0  <  .8  <  p   . 

The  well-known  relations  which  connect  the  components  of  a  vector  in  a 
paraboloid  system  of  coordinates  with  its  components  in  a  Cartesian  coordinate 
system  allow  us  to  obtain  from  (l,l)  boundary  conditions  for  the  required 
functions  u,v,   and  w.     Indeed,  since 


a  E     +  p(E  cos  0  +  E  sin  fi) 
z  X  y 


E0  =  E  cos  0  -  E^sin  0,     whereas    S     = 


relations  connect  u,  v,  and  w: 


"•"  -|P=Uo  ■  "le-P,- "iB-p,  *  ^p^^lP-P,  ■  °- 


We  represent  the  required  field  as  the  sum  of  the  field  of  the  source 
in  free  space  and  the  secondary  field, which  we  shall  seek  in  the  form  of  an 
integral  of  the  appropriate  sets  of  special  solutions  of  (1.2).  Thus  the 
complete  field  will  be  written  in  the  form 

E  «  E°  +  E-^,  E  =  E°  +  E""-,  E  =  E°  +  E"'-  , 
xxx'yyy'zzz' 

where,  in  the  notations  of  (l.l),  the  primary  field  (i.e.,  the  dipole  field)  is 
given  by  means  of  the  formulas: 


-  u  - 


"°-lfe 


-ikr 


£-5  (3  ♦  3ikr  -  k  r  )  -  (1  +  ikr 
2r 


-kV)], 


(1.5) 


-ikr 
op     e pz 


2  2, 


^Z_^-(3.3iicr-k'r'), 


^     ^-ikr      2  _  _ 

r  dT 

^  r  =  i  (a^  +  p2)^   p  oafi). 

In  the  case  under  discussion,  a  convenient  set  of  special  solutions  of 
(1.2)  regular  in  the  interior  of  the  paraboloid  is  the  set  of  solutions  of  the 
form 

\  -  ^(-^8^^^  (p)  a(^^)   (a). 


(1.6)  v^  -  a^(v)B^^^  (p)  A^^^  (a), 


where  A       (a)  and  B       (6)  are  the  wave  paraboloid  functions  connected  with  the 
hypergeometric  function  by  means  of  the  relations 


a[]Im   -  (ika2)2  e     "^F(v,.  .  1,   iKu^)  .  ii^'r^^,^^^^^^/^,..^^'^ ■ 
B^^jj|(p)  -  (ikp^)^e     ~^F(in+l-v,m+l,ikp^)  =  (ikp^)"^/\_(^^^)/2,m/2^^^^ 


and  satisfying,  respectively,  the  equations 


-  ^ 


A"*ir*  ^-2-2   ^^  •  .-^*1 


(1.6) 


B"^  1  B'*  [k2p2  .  2^  .  Uik  (H^  -  v]  ]b  -  0. 


Then  the  secondary  field  constructed  out  of  the  special  solutions  (].6) 
is  of  the  form 


u^  =  I  a^(v)  B,\-^;;  (p)  a),-'^   (a)  dv        0  <  R(v)  < 


J^^^  (p)  A^^) 
,T   o^  '   v.o  '■'^  ^     V,0 

^o 


(1.9)       v^  =  J  a^(v)  bJ,^^  O)  A^^^  (a)  dv       |  <  R(v)  <  1 


'  =  L^2^^)^i'2^-^)^l'2(»^ 


w  =  I   a^U;  b|^^^  (H)  A^^^  (o)  dv        1<  R(v)  <| 


where  the  paths  L  ,  L, ,  L_  are  parallel  to  the  imaginary  axis  in  the  plane  of 
the  complex  variable  v  (Fig,  2), 


t 

1 

^            1 

1 

\               / 

0 

1/2 

1 

3/^ 

Lo 

L, 

Le 

Figure  2 
The  choice  of  paths  of  integration  is  conditioned  by  the  fact  that  the 
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special  solutions  of  the  type  (1.6)  must  satisfy  definite  conditions  for  a  ->oo, 
nanely,  they  must  decrease  or,  at  least,  not  increase  as  a  ->oo  and  must  correspond 
to  a  wave  going  out  to  infinity.   The  unknown  functions  a. (v)  must  be  found  from 
the  two  boundary  conditions  (l.'4)  and  the  condition  div  E-*   0. 

2.  Determination  of  the  coefficients  of  the  representations 

In  order  to  find  the  values  of  the  coefficients  a.(v)  in  the  representa- 
tions of  the  components  of  the  secondary  field  it  is  necessary  to  make  certain 
assumptions  concerning  the  fonn  of  the  required  functions.   The  correctness  of 
these  assumptions  can  be  tested  after  we  shall  have  obtained  explicit  expres- 
sions for  those  functions.  Let  us  put 

(2.1)       a^(v)  -  -  ^  e-^V  (v)r  d  -  v)  c^  (v), 
a^(v)  .  ^  e-^"^  r  (v)r  (2  -  v)  c^(v), 
^^"^   "f^^""^"'  r  (v)r  (3-v)  C2(v), 

where  we  suppose  the  functions  c.(v)   regular  in  the  appropriate  strips  of  the 

plane  of  the  complex  variable  v  ■  er  +  ±T .     Thus,   c  (v)   is  supposed  regular 

(at  least)   in  the  strip  0  <  R(v)  <  1,   c-(v)   (at  least)   in  the  strip  -1  <  R(v)  <  2, 

and  c-(v)    (at  least)   in  the  strip  0  <  R(v)  <  2.     In  addition  we  assume  tliat  in 

these  strips  the  coefficients  c,(v)  are  of  the  order    Oit     e'  '^^o)  for 

T ->oo  and  of  the  order  ^i\T^\   for  f->  -oo. 

We  now  turn  to  the  condition  div  B^  0,     If  we  take  into  consideration 

One  should  keep  in  mind  the  asymptotic  behavior  of  the  functions  ^,      (a)  for 

5  v,m 

a  ->  oo:  V  1     °^     ito 

"'"^      a->oo   rCv)  '  P 

m  lioa'^ 

r(m+l-v)  ^  -^ 
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the  relatic.is  connecting  derivatives  in  the  paraboloid  system  of  coordinates 
and  the  Cartesian  system  of  coordinates,  namely 


8    cos  0  ff.    3 

■9  51    ^  =  3^^  ^  /^Q  9  4.  „  9^  4.  COS  0  a 


we  arrive,  on  the  basis  of  the  fomiilas  (1.9)   and  the  condition  div  E  -  0,  at 
the  following  equality. 


X%<v)foB<^)(„,i^^c,.^Bt^)-<^)0"' 


a  +p 


M:^<>^'S^^^-:i:^<A^<H 


(2.3) 


dv 


dv  + 


/    a,Cv) 


6  <:^'f>'v:r<") 


a  +p 


^Bi!r«"«i>>* 


^^« 


i^^^)  «ii^«'] 


dv  «  0. 


To  deduce  from  this  equation  a  connection  between  the  coefficients 
it  is  necessary  to  represent  the  integrals  in  (2,3)  as  integrals  in  the  same 
wave  paraboloid  functions.     Noting  that  the  integrels  under  consideration  are 
odd  functions  of  the  variable  a  we  use  for  this  purpose  the  functions  A^  -^(a). 
This  leads  to  the  following  integral  equation  for  the  required  coefficients  c^(v): 


c     (v)e-^''^  ,    c.(v)r-^''^ 


(2.U) 


_? , .  dv  -  2        4 , r  dv  *  2ie-^"^^  c,(n) 


If  we  write 


c.(v)e-^"^  r      c.(v)e-i''^ 


■inn  r    c,  ^v;e  y^      c  ^v;e 

2ie-^  c,(^)   -  j     4j^^rairvT  ^v  -   [       sin  n(^-vj 


dv  , 

then,  using  Cauchy's  theorem,  we  can  write  (2. u)  as  a  homogeneous  integral 
equation  with  the  kernel  cosec  n(vi-v).  This  leads  to  the  following  functional 
relation  for  the  coefficients: 

(2.5)     c^(v)  -  c^{v)   -  c^(v+l)  +  CgCv+l)  "  0. 

To  make  use  of  the  boundary  conditions  (l.li)  for  the  derivation  of 
other  connections  among  the  unknown  functions  c.(v)  it  is  necessary  to  represent 
the  primary  field  using  integrals  in  the  wave  paraboloid  functions.  Applying 
the  inversion  formula 

f(a)  -  —^^ f  e-^"^  r(v)r(m*l-v)f  A^^2(a)  dv, 

00 

f  -  J   <^f(^)A[^i  (a)  da      (^5  <  R(v)  <  2lij   , 


(f.6) 


"■it  is  also  assumed  that  n  is  located  to  the  right  of  L,  . 

A  formal  derivation  of  the  relations  (2,6)   is  based  on  1 

senting  the  integral  transform  under  consideration  as  the  result  of  repeated 

application  of  the  simpler  Mellin  and  Hankel  trfensfonns. 
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to  the  functions   (l.'j)  which  represent  the  field  of  the  dipolejWe  obtain  the 
representations 

"°  =  4  '°^'^'v!o^P>  ^i^o(^>  dv  (O  <  R(v)  <  i)  , 

o 

(2.7)  v°-J    aJ(v)B^^j(p)  A^^jCa)  dV  (l<R(v)<lj, 

w°  -  J    a°(v)Bj,^>(p)  A^^^(a)  dv  (l  <  R(v)  <  |)  , 


where 


°(v)  .  £^  e-^"^v2-v.l)  r(v)r^(i-v). 


o 

(2.8)  a°(v)  =  |g!  e-^"^(l-v)    pM    r^(2-v), 

-2(^>  =  -in«""'r(^^  r'(3-v), 

(3) 
and  where  B^"^  (p)  are  the  wave  paraboloid  functions  of  the  third  kind  connected 

with  Whittaker's  function  by  means  of  the  relation 

(2.9)  b(3)(P).(U*V^^W^_„,,^^(V). 

In  the  case  under  consideration  we  can  represent  both  sides  of  (l.U)  in 
the  saine  paraboloid  functions  using  the  recursive  formulas 


(2.1D) 


>^  "I'llf")  ■  (■»*!) [4'i(°>  -  *v-l,„(«)]' 
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for  the  functions  A^  (a),  which  leads  to  the  following  relations  among  the 
required  coefficients  J 

r^o  W^Mllo^h^  -  ^o^^^^v'o^Po)]  '^y^^^^'^'^^'Mll^^' 


.  c^iv.lMl.v)Bliy^^)]   -  -r(l-v)Bi3)   (p^)  -r(3-v)B2i   (p^). 


(2.11) 


=o(^);^o^l'o(Po)  -  ^l(-)(l-)Bi;j(p^)  .  c,(v-l)vBi^^^,(p^) 


-  r(2-v)vBj,3)(p^) .  ri2-v)Blil^{^j. 


Making  use  of  (2.?),  (2.11),  and  function  theoretical  considerations, we 
arrive  at  the  equation 


(2.12)    c^(v)  -  c^(v-l)  =  f(v). 


Tirtiere 


'^•»>     ''"' ,a)^„ m' 


o b^ 2 ^ ^ i i 

(2-v)Bi^i,l(pJ  -  (1-v)b(^^  (p^) 


We  note  that  the  function  f(v)   is  regular  at  least  for  R(v)  <  2,  since 
both  numerator  and  denominator  are  entire  functions  in  the  part  of  the  plane 
under  consideration.      In  addition  the  expression  in  the  denominator  does  not 
vanish.     If  we  were  to  assume  the  opposite,   then,   as  can  be  shown,  the  con^s- 
ponding  homogeneous  problem  would  have  a  non-trivial  solution;  a  conclusion 
which  contradicts  the  uniqueness  of  the  solution  for  the  class  of  problems  under 
consideration. 
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The  solution  of  equation  (2.12)  is  of  the  foim 

00 

(2. Ill)      c^(v)  =X  f(v-n). 


Using  (2,5)  and  (2.11)  we  obtain 

,  ,     ,  ,   l-vr(l-v)[(l-v)B'^j(P„)BW(p  )  ■  B^|)   (p^)B^^>   (p^)] 

C^(v)  =  C^(v)  +  2 2 i^=2 2 2 

(1-v)b(^>  (Po>  *  ^^111,1%^ 


(2.1$) 


l.(l-v)r(2-v)[(2-v)B21^l(P,)B^^j^,(p,)  -  B(3)(p^)Bg(p^)] 
(v) !:— 75 ^ Hj * " 


C2(v)  «  c^^,v; ^ 2" 

3.  Investigation  of  the  solution.  Transition  to  geometrical  optics. 

In  view  of  the  results  of  Section  2,  the  final  expressions  for  the 
components  of  the  secondary  field  are: 

^x  -=  -  Si/,  %(-)«■""  r(v)r(l-v)B(l>)  Alalia)   dV 
o 

*^/  C2(v)e-^"T(v)r(3-v)Bi^^(?)  /^^(a)  dv  cos  21^, 

^2 
(3.1) 

^y  "  ^  /j^  C2(v)e--T(v)r(3-v)Bi^^(p)  ^^^^(a)  dv  sin  20, 
e1  =  1^1  c^(v)e-^"T(v)r(2-v)Bi^^(p)  /^^^(a)  dv  cos  0, 

*It  should  be  noted  that  the  solution  of  (2.12)  is  defined  uniquely  to  within  a 
periodic  function  with  period  one.  On  the  basis  of  the  assumptions  about  the 
coefficient  c, (v)  it  can  be  shown  that  in  the  case  under  consideration  this 

periodic  function  vanishes  identically. 
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where  the  coefficients  c^^Cv)  are  defined  by  the  formulas  (2,13)-  (2,l5). 

If  we  assume  a  <  p,  then  it  is  possible  to  supplement  tiie  paths  in  the 
above  integrals  by  arcs  of  circles  of  radius  R  ->oo.  Inside  the  closed  contour 
the  integrands  are  mesomorphic  functions  with  poles  at  the  points  v  =  -  m 
(m  =  0,1,2,  ...),  Using  the  residue  theorem  we  find  the  com.ponents  of  the 
secondary  field  in  the  fonri  of  Pinney  series.  Thus,  for  example, 


(3.2)    eJ  .  kis! £  ('"^i^B^iJi^p^^-i,^^)^^^  <^  ^  ^^-p)' 


2  00  ,^y,  /,  ^  00 

z 

p=m 

where  the  functions  •'^l™  n('i)  are  connected  with  the  Laguerre  polynomials  by 
means  of  the  relations s 

ika^  .1,..    2x 

These  series  converge  only  in  a  part  of  the  domain  (for  c  <  p)  so  that 
one  cannct  vrithout  further  proof  claim  that  they  represent  the  required  solu- 
tion. 

Our  solution  of  the  problem  in  the  form  of  contour  integrals  (3.1)  can 
be  used  for  the  investigation  of  the  electromagnetic  field  at  any  interior  point 
of  the  paraboloid  in  many  special  cases.  In  this  paper  we  shall  dwell  on  the 
limiting  case  of  very  short  waves  when  ka  ->oo  (a  is  the  distance  from  the 
vertex  to  the  focus  of  the  paraboloid). 

In  view  of  the  asymptotic  representations  of  the  wave  paraboloid 
functions  for  large  values  of  the  argument,  the  coefficients  c.(v)  have  the 
following  form  as  k  ->oo: 


Jur  results  show  that  the  Pinney  series  represent  the  solution  in  that  part  of 
the  interior  of  the  paraboloid  in  whiich  they  converge. 
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e.(v)e--^   =     r(^-")r(l-v)  (kp^)2^  ,""'^° 


°  k->oo         ik^^ 


(3.U)  c,(v)e-i-  ;^     inO-jOt^  (,p2)2v 


•ikp? 


k->oo       (ikp    )- 


e  rv)e-i-  ^     r(3-v)r(2-v)   ^     2)2v  ^-^^0^ 

2  k->oo     (ik^;)3  ° 


utilizing  the  integral  representation  of  the  product  of  two  wave 


paraboloid  functions 

inm       k/2     ik/   2  q2 


v,m'      v,m  r(v)rU*l-v)  J^ 


X     «J  (kap  sin  2e)d9,  (0  <  R(v)  <  m+l), 


and  the  asymptotic  formulas  (3.  U)   for  the  coef ficients^ we  can  transform  the 
integrals   (3.1).     In  this  way  we  obtain  for,   say,  the  component  E     the 


relation; 


•ikp^      ^n/2         ik(2  ^,2^ 


e1     ^    .8£^«_  e      2  J  (kaB  gin  29)-^ 

z  .    .  cm     ,  „6  1  ^_2a 


^k->oo"^"i"^      i  -  tg-e 


X      f    (2-v)r^(2-v)(kp2^g  e)2^  dv  cos  0. 


In  view  of 
r    (2-v)P^(2-v)(kp^tg  9)^^  dv  =  UnikSo°tg%  K^(2kp^tg  8), 


Hi 


where  K^(x)  is  the  McDonald  function,  we  obtain,  upon  putting  tg&  -  x, 


,2       00       ikf2  „2vl-3C^ 


,}  ^  .32£UflU  -''''"or   ■-<''-'"ii?.y2taej\K,,«2.)i^ 


k-^-oo       ° 


K^pu  ^-^"  o  r  ,-  -^°  -"  't:^  ji^\  ,^(3,p2^,  ^ ,,, 


2 

If  we  introduce  a  new  variable  2kp  x  =  t  and  note  that  for  k  ->oo  the 

o 

principal  part  of  the  integral  corresponds  to  values  of  x  close  to  zero,  then  we 


obtain,  approximately, 
—(a 
^i  -  ''-^'^ e-^°co3^|   Vt).i^ 


.^(aV)   ..o2 


or,  noting  that  (cf.  [3]  for  example) 

00 

r  K  fatU  fh+^t^^*^-'^dt  =  (2a)^(2b)^  P(n^v^l) 


ik,  2  o2-   -ikpf 
E^  ^  .l^H^-r^'^-^  '  e    °  gp    cos 


c  °         ,u   „2  r 

k->oo  '^     " 


1  +, 


3  • 


Proceeding  in  an  analogous  manner  in  the  case  of  the  remaining  components 
of  the  secondary  field, we  obtain,  under  the  assumption  that  ka  ->oo,  the  follow- 
ing asymptotic  formulas: 
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Kf]' 


-ik(z+2a)       TT-  cos  0 


1      ^    _  16£  e       ^  2a 


k->oo       °  (2ar 


T     2     2  2 

where  we  have  put  z  =  ^(a  -p),   p=ap,p     =2a. 

These  approximate  formulas  for  the  components  of  the  secondary  field 

admit  of  a  simple  physical  interpretation,   since  they  represent  the  field  of  a 

non-homogeneous  plane  wave  whose  phase  and  amplitude  satisfy,   as  is  easily  seen, 

the  requirements  of  geome+rical  optics.     Indeed,   the  formulas  (3.5)  show  that 

as  k->oo  the  z-component  is  of  a  higher  orxier  of  sma3.1ness  than  the  components 

It    and  ir  .     If  we  suppose  the  electrical  field  known  and  compute  the  magnetic 
X  y 

field  by  means  of  the  formula  H  =  -  -«-  rot  E,  we  arrive  at  the  conclusion  that 

here  too,   as  k->oo,  there  remain  only  the  x-  and  y-components  of  the  field ^ -^nd 

that  2 

^-ik(z+2a)  [^)     sin  20 


(3.6) 


„1  2pk  e 

„1     _        2PIC  e-^'^*^^'     l-(fef  cos»i 
Thus,   in  view  of  the  formulas  (3-5)  and  (3.6),  the  mean  value  over  a 
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period  of  the  only  non-trivial  component  of  the  Poynting  vector,   namely,   the 
z-component,   is  given  by 


(3.7) 


dS 


pdpd0. 


On  the  other  hand,  if  we  compute  the  stream  of  energy  radiated  by  the 
dipols  and  contained  in  the  cone  between  the  angles  9  and  9+d9  and  if  we  assume 
that,  in  accordance  with  the  laws  of  geometrical  optics,  this  whole  stream  raustj 
after  reflection  from  the  paraboloid,  find  itself  in  the  cylindrical  layer  of 
thickness  dp  (Fig.  3),  where  dp  is  connected  with  the  angle  d©  by  means  of  the 
relation 

ri©  =  ^ . 


[Hm 


we  are  led  to  the  equality 

2  2 
(3,8)  dw  =  £-^  sin^  sin  e  d9d5^, 

OTTC 

where  Y  is  the  angle  between  r  and  the  x-axis,  i.e.,  cos  T  =  sin  9   cos  0. 


i 

^x 

;-^ 

( 

^' 

t 

\ 

\ 

\ 

0 

z 
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Fig\ire  4 


Noting  that  sin  © 


{-(^. 


we  obtain  the  formula 


(3.9) 


dw 


8nc?  T77p¥r 


^^^i^ 


which  coincides  with  the  equality  (3»7)  obtained  directly  from  the  asymptotic 
expressions  for  the  electromagnetic  field  inside  the  paraboloid  for  k->oo. 

The  fact  that  the  phase  of  the  secondary  field  is  also  in  agreement  with 
the  laws  of  geom.etrical  optics  is  easy  to  verify.  Indeed,  if  we  consider  in 
Fig.  h   the  paths  of  rays  emanating  from  the  radiating  dipole  located  at  the 
focus  of  the  paraboloid,  then,  according  to  the  laws  of  geometrical  optics,  the 
phase  of  the  secondary  field  at  the  point  M  will  be  k(0/.  +  M).      T.-^  we  use  the 
geometric  relations  which  obtain, we  find  that  OA  =  AB  =»  2a  +  z. ,  so  that 
OA  +  AM  »  2a  +  sl.,  i.e.,  the  seo<->ndary  field  phase  is  exactly  equal  to  the  phase 
which  enters  into  the  relations  (3.5)  artd  (3.6), 

In  conclusion  the  author  wishes  to  thank  G.  A.  Grinberg  for  having 
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suggested  the  topic  of  this  paper  and  to  N.  N.  Lebedev  for  many  helpful 
suggestions. 
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APPENDIX 

W,  Magnus  and  N.N,  Lebedev 


Since  it  is  not  obvious  how  to  derive  equation  (2.L)  in  this  paper,  a 
proof  of  this  equation  will  be  sketched  here.  We  are  indebted  to  Professor 
N,  N.  Lebedev  of  the  Physico  Technical  Institute  in  Leningrad  for  communicat- 
ing this  proof  to  us  after  consvdting  with  the  author  of  the  paper. 

The  derivation  of  (2,U)  is  based  on  the  assumption  that  a  solution  of 
the  equation 

regular  within  the  paraboloid  and  satisfying  the  radiation  condition 

(A2)   V  -  0  (a'*),  E  ^  *  ^^  '  °^°-''^^  (a->oo) 

can  be  expressed  in  the  foim  of  an  integral  of  the  type 

(A3)   V  -  v(a,p)  =  j  MCti.)  bJ^^|(P)  A^^^(a)  dji, 

h 
In  order  to  determine  the  coefficient  M(n)  let  us  divide  (A3)  by  /He  { 

and  pass  to  the  limit  p->  0. 

This  gives 

(All)    f(a)  -  llm  ^^5^  -  f  M(tx)  t^l^ha)   dji, 
p->  0  /lie  p    -'l 

from  which  it  follows,  in  consequence  of  the  inversion  formula  (2,6),  that 

00 

(,5)   M(^)  -  ke---^  ^—j-   f^'^)  i]l^-^  ^-' 
o 

It  may  be  easily  proved  that  the  functions  in  brackets  under  the 
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integral  signs  in  (2.3),  and  which  will  be  denoted  in  what  follows  by  v  (a,p), 
'^l(°-f?)t   V2(a,p),  are  regular  solutions  of  equation  (l).  v  (a,p)  and  VpCijP) 


also  satisfy  condition  (A2). 

Usi  ng  this  method  we  get 


(A6)   v^(a,p)  .  ^  b(1>(p)  ^]l\a)  .  -^  B^l)'(p)  A^^^Ca) 

f^ia)   -  li:„  ^l!ii.  -  -^  Ul^l'i^)  -   (l-2v)ika  ^l^lia)] 
°    p->  0  /^  p     yn?  a^  L  v,o  v,o   J 


.  v(l-v)  2^  ^illoM. 


from  which 


(A?)  M^(n) 


^ -(i-''-'""  iS^  f '^U'-' *i>)^- 


Taking  into  consideration  that  0  <  R(v)  <  »  and  ^  <  R(n)  <  1  and  making 
use  of  the  integral  representations 

A(l)(a)  =  2ei"A  ^^   e'^'/^  ^^^^^   at)e-^'t2.-2, 

^'^  r(n)  i^  ^ 

(A8) 

we  obtain  after  same  calculations  which  we  shall  give  at  the  end  of  our  analysis: 
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Therefore 

(Alo)     M  (X)  -  k  4=ii rSin-nv 

o^    '  sin  nti     sin  n(,|j.-vj 

The  case  of  the  function  v^Cajp)  may  be  treated  in  a  similar  way. 
The  function  T,(o,p)  does  not  satisfy  condition  {1^2).     But  if  it  is 
written  in  the  form 

(All)   v-^(a,p)  =  y*(a,p)  -  iJcB^^jCp)  A^^^(a) 

the  method  described  above  may  be  applied  to  v^(a,p)  where  one  must  put 

I  <  R(v)  <  R(n)  <  1. 

Introducing  into  (2.3)  the  expressions  obtained  for  the  functions  v  (o,p), 

T-(a,p),  v«(a,p)  and  equating  to  aero  the  coefficient  of  B^   l{f>)  tr  iMt  we  get 
■i-  c.  ^jJ-  !-'■»-'• 

the  integral  equation  (2,h)» 

We  still  have  to  derive  equation  (A9)  from  equations  (A8).     This  can  be 
done  in  the  following  manner:     We  have  from  (A8) 

00 


I    A::^,^   ,(a)   A:,%'(a)dfl 


o  -inA     in(ti+v)/2      /°°       ^^  r"^  .,^2   .2^  „„   ^  ,    ^^ 

^®     -2 da      J,(2Aat)J,(2Aat)e^^^  -*  V^'~'■t^■^''dtdtr 

r(2-v)        J   ) J„  ^     2 


r(^^).  -    -      -o     o'o 

We  shaU  exchange  the  order  of  integration,  carrying  out  the  integration  with  respect 
to  a  first.  The  well  known  Sonine-Schafheitlin  formulas  show  that 

,00  r .  ^'2 


r  1  t  'C-'^  if  0  <  t  <  r 

J  J^(2^at)  Jo(2)1^aT)do 


0    otheivise. 


-  22 


ITierefore,  we  find  for  the  left-hand  side  of  (A8)   the  expression 


:l^l^^^    far  L  ,ia^-t^),3.-H-i- 

(2-v)2K1c         /  J 


,,,,,        8g-inA-in(t..v)/2      ^  f        i(^2_^2)   ^^^     ^^_2V 


r(n)r(2-v)2)^ 


2 

By  expanding  e~^       in  a  power  series  and  integrating  term  by  term  with  respect 

to  t  we  find 

r 

(A13)  1   e-^*^  t^^^-^  dt  =  4^  ,F,    (nj  ^*l;  -it^), 


;a13)        j^e-^*^  t^^^-^  dt  =  -|^  ^F^  (nj  ^*l;  -it' 


o 
The  second  integration  in  (A12)   reduces  now  to  an  evaluation  of  the  integral 
2 


(AlU) 


l_j    e^-^V^^-^^-^F,  (.;..l,  .ir2)dr. 


By  changing  the  path  of  integration  or,  equivalently,  by  substituting  a  new 
variable  s  for  -iT  ,  we  find  for  (Alii); 


exp[in(^.-v)/2]|  ^-s  ^^x-v-l  ^^^^^,   ^^^,   ^^^^ 


^    r(^.-v)exp[in(^.-v)/2l^^_^  (^.^^^.  ^,,.  ,) 

_.  (i^)-l  r(^^-v)r(^i^l)r(v^l-|x)   exp[in(^..v)/2] 

rci)  rcv-^i) 

o  "  ^(^t)    exp[iTi(^i-v)/23 
H  ^^^^^j       3in  it^n-v;   • 

Combining  this  result  with  (/13)  and  (A12)  we  find  that  (A9)  is  true. 
The  analysis  given  is  purely  formal. 
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Foreword 

The  following  paper  by  Kontorovitch  and  Muravev  is  concerned  primarily 
with  the  asymptotic  evaluation  of  double  integrals  of  the  form. 

S 

Because  this  subject  was  being  investigated  by  two  of  our  own  people  and  it  was 
desirable  that  we  learn  the  contents  of  the  Russian  paper,   a  translation  was 
undertaken. 

The  purpose  of  the  paper  is  to  extend  to  double  integrals  the  application 
of  the  method  of  stationary  phase  which  in  the  past  had  been  used  to  evaluate 
single  integrals.     The  authors  do  make  such  an  extension  in  that  they  obtain 
the  principal  part  of  the  double  integral  and  show  that  it  is  of  the  order  of  l/k. 
Their  result  is  shown  in  formulas  (9)  and  (11)  of  their  paper.     Actually  the 
term  obtained  by  the  authors  is  only  one  of  a  number  of  terms,   some  of  the  same 
order  and  others  of  higher  order.     These  other  terms  come  from  other  types  of 
stationary  points,  or  critical  points  as  they  are  called  in  the  literature,  but 
these  are  not  considered  by  the  present  authors.     They  limit  themselves  specifi- 
cally to  the  condition  that  the  phase  function  9^  contains  at  most  one  stationary 
point  and  moreover  a  stationary  point  at  which  the  phase  function  is  either  a 
maximum  or  a  minimum. 

The  results  obtained  by  these  authors  have  now  been  superseded  by  the  work 
of  a  number  of  others.     I  should  like  to  refer  to  our  own  Research  Report  No.  EM-100 
by  Douglas  S.   Jones  and  the  undersigned  entitled   '/Asymptotic  Exransion  of  Multiple 
Integrals  and  the  Method  of  Stationary  Phase.'     This  report,  as  well  as  a  number 
of  others  cited  in  the  introduction  and  bibliography  indicate  the  more  recent 
work  that  has  been  done  and  the  proq;re3S  made  beyond  that  in  the  paper  by 
Kontorovitch  and  Muravev. 

Morris  KlirB 
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Abstract 
In  this  paper  we  consider  the  deirivation  of  an  asymptotic  formula  for 
surface  integrals  of  oscillating  functions  (the  method  of  stationary  phase). 
3y  means  of  the  fonriula  so  obtained  we  establish  a  connection  between  the  laws 
of  reflection  of  geometrical  optics  and  the  laws  of  classicrl  electrodynamics 
taking  as  our  starting  point  the  integral  formulation  of  the  diffraction 
problem. 

1.  Preliminary  remarks 

h   large  number  of  problems  in  mathematical  physics  and,  in  particular, 
in  radio  physics  and  in  optics  involves  the  computation  of  integrals  of  fast 
oscillating  functions.  The  study  of  such  integrals  is  of  interest  not  only 
from  the  point  of  view  of  solving  concrete  problems  but  also  from  the  point 
of  view  of  establishing  general  relations  which  are  valid  in  the  case  of  very 
high  frequencies. 

In  particular,  it  is  possible  to  show  that  the  problem  of  connecting 
many  of  the  laws  of  geometrical  optics  with  classical  electrodynamics  can  be 
reduced  to  the  problem  of  studying  the  behavior  of  integrals  of  the  type 
mentioned  above  in  which  the  integrands  are  functions  with  large  frequency  of 
oscillation. 

For  one -dimensional  integrals  of  this  type  the  corresponding  asymptotic 
representations  can  be  obtained  using  a  well-known  formula  which  is  a  consequence 
of  the  principle  of  stationary  phase. 

However,  it  is  not  possible  to  justify  fully  the  extension  of  this 
formula  to  the  case  of  integrals  taken  over  a  two-dimensional  domain  and  it  is 
precisely  domains  of  this  type  that  usually  turn  up  in  the  problems  mentioned 
above. 

With  these  considerations  in  mind  we  give  below  the  derivation  of  an 
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asymptotic  formula  which  permits  us  to  obtain  approximate  expressions  for 
integrals  of  the  type  mentioned  above.  This  formula  represents  an  extension 
of  the  method  of  stationary  phase  to  the  case  of  two-dimensional  integrals. 

Since  a  surface  integral  can  be  reduced  to  one  or  more  integrals  over 
plain  domains,  we  restrict  ourselves  to  the  derivation  of  an  asymptotic  formula 
for  a  plane  dornain.  The  result  obtained  can  be  extended  to  the  case  of  a 
curved  surface. 


2.  Derivation  of  the  asymptotic  foimula 
Consider  the  integral 


(1)       1=1  Fe^  dS 
S 


/.' 


taken  over  a  plain  domain  S.  Here  F  and  0  are  real  functions  of  two  variables 
X  and  y  defined  in  the  domain  S  and  on  its  boundary  P  and  k  is  a  constant  real 
number  which  we  shall  assume  to  be  large. 

We  impose  the  following  restrictions  on  the  functions  F  and  ^:  F  is 
single-valued,  continuous  and  finite  together  with  all  its  derivatives  of  first 
and  second  order.  ^  is  single -valued,  continuous  and  finite  together  with  its 
partial  derivatives  up  to  the  fourth  order,  and  the  interior  of  S  contains  at 
most  one  stationary  point  in  which  grad  {(T  '  vanishes  and 


a^czf  a^rf    (^ 


(2)     f^  ^  -  life)  >  ^' 


Obviously  sufficient  but  not  necessary. 
**In  the  sequel  grad  and  div  denote  surface  differential  operators, 


i.3,,   the  stationary  point  corresponds  to  a  minimum  or  maximum  of  the  function  0. 


Let  Y  be  the  boundary  of  S  (cf.  Fig.   l).     We  have 


Jre^dS     +    I    Fe^ 


ik0 


It  is  easy  to  obtain  the  identity 


(3) 


-^--i 


e^  div  -JSl 


•^     il(0 


div 


F  q  e 


(q  grad  0)  (q  grad  0) 


where  q  is  an  arbitrary  vector. 

/.t  fii^t  we  assume  that  nowhere  in  the  domain  under  consideration, 

except  for  the  stationary  point,   is  grad  0  perpendicular  to   the  radius  vector  r 

connecting  the  stationary  point  with  a  point  of  S. 

Then,  putting  q  =  r  ,  where  r    ■  ,  we  can  write: 

^  ^        0'  o       ,-, 


I  =   f  Fe^  dS  -  ^         f    e^  div  J^ dS  -   f      di 


Fe^r. 


dS 


sis  (Vrad^)  ^_5         (r^grad0)      _ 

o  o 

or,   on  the  basis  of  Ostrogradsky's  theorem. 


(U) 


-  fpe^dS  -i.  [     r   e^div^'^°  dS-   r   I^!^(iA)  dl  - 

/g  "^   [  sis  (^^grad  />)  J^  (r^grad  0 ) 


■{ 


Fe 


^(1^^ 


dl 


Y  (r^grad  0) 


Here  m  and  m     are  unit  vectors  of  the  normal  to  the     contours    |      and  y 
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Now,  let  the  radius  r  of  the  circle  y  tend  to  zero.  In  view  of  the 


^   j,      (^  .rad  0)     J 


.^ft 


2_  dl  +  lixn    -^TT—  dl 


r  (V^d0)       r^O.^  ^ 


Using  polar  coordinates  r  and  9  with  origin  at  the  stationary  point  and 
ignoring  quantities  of  order  higher  than  r  we  can  write 


^  =     a-cos^  +  2a_cos  6  sin  9  +  a  sin  9  r 


where 


'■i  .r--^  .:•'■& 


st 


Thus 


lim  I|^  dl  =  F(x  .y  )e  °     °  . 5 


de 


5- 

pcos  &  sin  9+a-sin  9 


Under  the  assumption  (2),  i.e»,  for  a.a  -  a.  >  0,  this  integral  can  be 
computed  by  means  of  the  substitution  tg©  =  z  and  is  equal  to 


I 


2n 


d9 


a,cos  9+2a2Cos  9  sin  9+a-.sin 


5" 


/". 


dz 


2n 


-00   123    /  i^'i  ~  °2 


We  thus  obtain  the  following  expression  for  I 


5  - 


(5)      I  =  - 


H 


2n  ,,  .^^V^^o^        ^Fe^^U^^) 


(r^grad  0) 


dl  + 


I 


.       e^ 


Fr. 


(r"  grad  0) 


dS 


This  implies 


(6) 


ik  I 


2n  .,       ^^Wt'^o-yo)       f   fe«'(»?,) 


^ 


:^  F(^o>yo)« 


fe^div- 


,grad  i^) 


dl 


r 


^. 


,grad  0 ) 


We  now  consider  the  integral  on  the  right  side  of  (6) 


^   =  fe^  div    ^    '^°         dS  .     f   e^  div  - 
{  (r.grad^)  3^^  (^^ 


Fr, 


grad  ^) 


dS  =  I     +  I 


where  S^  is  a  circular  domain  of  small  radius  r,   about  x  ,y  . 
o  00 

Writing  the  integral  I     in  the  form 
2n    .r. 


/^^ 


-   e^ 


Fr. 


I  -  I  e"""^  div 

S 


(r^grad  0 ) 


dS 


o    0         \ar  / 


dr  d9. 


we  shall  show  that  the  integrand  in  I  is  finite  in  the  whole  domain  of 
integration.  He  can  write 


9? 

(f)- 

For  small 

t 

r 

Br 

'^ 

1^  =  (a^cos^e  +  2a2Cos  6  sin  9  +  a  sin^9)r  +  O(r^) 

so  that,  in  accordance  with  the  inequality  (2),  the  first  term  remains  finite 
for  all  e. 

As  for  the  second  term,  we  can  write 

^^    6r^    dr^   ^  8r^    br^  Br^  ^  ar^ 


It  is  now  easy  to  see  that  ^  ["^t)   ^^   bounded  in  all  of  the  domain  S. 
We  can  therefore  write       V 


Is    de 


_a_   Fr 

3r  "^ 

3r 


dr. 


If  we  now  denote  by  M  the  largest  value  of  the  integrand  in  the  domain 
of  integration,  then 


(7) 


I  1  <  2nr^M 


and  M  is  independent  of  k  and  r, . 

We  now  transform  I  using  identity  (3)  and  replacing  F  in  (3)  by  F^j 


7  - 


Fr. 


(r^grad  0) 


Then 


ik0    ,.          ^1^0  ,„ 

e    ^  div  — r dS 


/S 


F^e^(J?r^) 


dl 


S-S, 


(r^grad  0)  r->       (r^grad  0) 


Since  the  integrand  remains  finite  in  the  domain  and  along  the  paths 
of  integration,  we  may  assume  that 


(8) 


li-Kf 


where  N  is  a  constant  which  is  independent  of  k  but  dependent  on  r, . 

Given  c  >  0  we  can  make  |l  |  <  e/2  by  choosing  r,  small  enough. 

Now,  keeping  r^    fixed  in  accordance  with  (8)  and  choosing  k  large 
enough  we  can  make  ll  |  <  e/2. 

Thus,  for  a  sufficiently  large  k,  the  right  side  of  (6)  can  be  made 
arbitrarily  small. 

Hence  the  formula: 


(9)   (pe^dS 


1 
Ik 


/ 


Fe"^rtr.) 


dl 


^hO(h. 


the  expression  (r  grad  0)  does  not  vanish  in  the  domain  of  integration.  We  shall 


now  try  to  remove  this  restriction. 

In  view  of  the  inequality  (2)  we  caai  always  find  in  the  vicinity  of 
points  X  ,y  a  finite  circular  domain  <J"  with  boundary  r^   in  v/hich  (r^grad  0)  /  0. 
Then,  using  (9),  we  can  write 


ikJiJ 
^ F(x  .y  )e    °  °  -   ~ 


2n    „,_  _  .."-f^f^.^o'   f   "'^'Vo)^ 


where  m,  is  the  outer  normal  to  y-i  • 

We  now  consider  the  remaining  part  of  the  domain  of  integration  S-S 
bounded  by  the  contours  P  and  f-,   and  not  containing  stationary  points.  Next 
we  choose  a  vector  q  (which  appears  in  formula  (3))  in  such  a  way  that  q  is  not 
perpendicular  to  grad  0   on  Yt  and  is  equal  to  r  on  f-,  • 

Before  we  continue  our  argument  we  show  that  such  a  vector  can  at  all  be 
constructed. 

We  choose  q  »  F^r^  +  Fpgrad  0,  where  F^  and  F^  ht^   functions  of  x  and  y 
such  that  F^  =  1  and  Fg  -  0  on  y^.     Then  (q  grad  0)  "  F^  H  +  FgCgrad  j6)  and 
q  «  r  on  Yt  In  view  of  (2),  |=-  does  not  change  sign  inside  c  .  Hence  the 
nearest  point  at  which  ^  =  0  lies  outside  the  contour  Yi • 

We  can  thus  select  a  domain  cr  ^  with  boundary  Yo  containing  the  domain  cr 

such  that  ^  does  not  change  sign  in  o-,  .  We  assiune  that  F„  is  positive  every- 
or  -L  t 

where  (except  for  the  contour  Yt )  and  that  F,  has  the  same  sign  as  ~  has  in  the 

*We  assume  that  q  satisfies  the  sane  conditions  which  the  conditions  imposed 
on  0  imply  for  grad  0. 
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vicinity  of  the  stationary  point. 


(q  grad  jf)  -  F2(grad  0)^  +  ^^  ^  >  0- 
In  the  domain  exterior  to  Yo  we  have 

KtgradgT)!   >  FgCgrad  0)2  -   |f^  |^|    . 

Since  in  the  domain  under  consideration  F^  and  grad  0  do  not  vanish, 
it  is  always  possible  to  choose  F^  so  that  F^Cgrad  0)  >     F,  ^   .       Consequently, 
(q  grad  0)  >  0,  which  is  what  we  set  out  to  prove. 

Now,  using  (3)   and  the  theorem  of  Ostrogradski  we  can  write: 


/ 


Fe^dS  =  - 


H 


fe^div-Ia—dS-f       l£^^ll 
f_^  (q  grad  0)  p^     (q  grad  0) 


dl 


/.s  above,   applying  formula  (3)  once  more,  we  obtain 

ds  =  Oih. 


j  e^div  _la 


S-o- 


(q  grad  0) 


We  thus  obtain  the  following  representation  for  the  integral  in 
question: 

One  way  of  doing  this  is  to  put  F_  =  M(r-r  ),  r  =  distance  from  the  stationary 
point,  M  -  a  sufficiently  large  positive  number,  r  =  radius  of  the  circle  f-,, 
and  to  put  F  «  1. 


ID  - 


(10) 


Jre^dS-I,.!,.-^ 


/v 


^^(V.o)^"''^"*''"' 


,^ft  ^^ 


r  Fe^(mt) 
;  (q  grad  0) 


dl 


3^2 


^IO(^). 


Depending  on  the  definition  of  q  on  the  boundary  of  the  domain  we 
can  obtain  different  representations  of  the  integral  under  consideration. 

In  particular,  if  we  put  q  =  grad  0  in  the  contour  integral,  then  we 
find: 


(11)   J 
S 


Fe^dS 


1 
IE 


/o^Oj^l  J     (grad  ^r 


*io<i). 


If  we  assmne  that  q  =  r  on  |  ,  then  we  again  obtain  formula  (9). 

It  is  easy  to  see  that  if  9^  is  constant  on  the  whole  contour  or  on 
each  distinct  part  of  the  path,  then  the  contour  integral  in  the  formulas  (9), 
(10),  and  (11)  is  of  the  same  order  with  respect  to  k  as  the  first  suiranand. 

In  the  remaining  cases  this  integral  will,  in  general,  -tend  to  zero 
with  increasing  k. 

It  is  easy  to  extend  formula  (9)  to  the  case  when  the  surface  is  curved 
and  is  characterized  by  Gaussian  coordinates  u,  v,  and  when  we  are  given  the 
values  of  the  second  derivatives  of  the  function  0  with  respect  to  these 
coordinates  at  the  stationary  point: 


n 


5^ 


St 


St 


Let  o-  be  a  plane  domain  which  is  the  projection  of  S  onto  a  certain 
plane  and  let  \       be  the  projection  onto  the  sane  plane  of  the  boundai7  P  of 
S,  Then 


:i2)   f  Fe^l'dS  -  f  Fe^H(ic,y)do-  ■  -  ^ 


2l«£^yU,y)e"*"'"'"°' 


Ft 


^A 


-/_ 


>ik0^:^ 


FH(x,y)e^(m^r^) 


dl 


(r^grad  0) 


*\0^\), 


where  m  is  the  nonnal  to  \      .   r,  is  the  radius  vector  on  cr  with  center  at 
the  projection  of  the  stationary  point,  and  H(x,y)  is  the  discriminant  of  the 
first  quadratic  form  of  the  surface  in  the  x,y  coordinates. 
Noting  that 


fTA  ■  f^^^  \^  ' 


where  ^r^^    {       is  the  Jacobian  of  the  coordinate  transformation  from  x,y 
to  u,v,  we  find  that  the  quantity 


H(x,y)  _  H(u,v) 


where  H(u,v)   is  the  discriminant  of  the  first  quadratic  form  of  the  surface  in 
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the  u,v  coordinates,  is  invariant  with  respect  to  any  coordinate  transformation. 

The  connection  between  an  element  of  arc  length  dl  of  the  contour  P 

o  '  o 

and  an  element  of  arc  length  dl  of  the  contour  \     is  given  by  the  relation: 


dl. 


(Vo) 


dl. 


where  m  is  normal  to  the  contour  |  . 

Thus  the  representation  (12)  can  be  written  down  as  follows: 


/■ 


(13)  1  Fe^dS  -  -  ^ 


2nH(u^,v^)        ^(^o^^'o^ 


/ 


a^B^-a^ 


-/ 


Fe 


^(m? 


(r^grad  0) 


dl 


*^0(i). 


If  instead  of  the  function  F  there  is  given  a  vector  function  F,  then, 
applying  (13)  to  each  of  the  scalar  components  of  the  integral,  we  obtain  the 
following  representation  for  an  integral  of  the  type  under  considei^tion: 


(U.)   JtelWdS  .  -  ^ 


f 


-J  0      0 

^^3-^2 


^iO(|). 


.     fe^(^r^) 


dl 


Here  the  tenti  r  (l)(r)  stands  for  the  expression 


(r  erad  5^) 


FdS. 
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The  above  method  for  obtaining  an  asymptotic  representation  for  surface 
iiitegials    [forrrulas   (9)-(lii)J    is  in  fact  a  generalization  of  the  so-called 
method  of  stationary  phase  which  is  widely  applied  for  approximate  computation 
of  one-dimensional  integrals  of  oscillating  functions. 

3,     Derivation  of  asymptotic  formulas  for  the  computation  of  the  electromagnetic 
field  reflected  from  a  surface  of  arbitrary  shape 

The  above  fon^oilas  can  be  used,   in  particular,   to  establish  a  connec- 
tion between  the  laws  of  geometrical  optics  and  the  laws  of  electrodynamics  and 
to  obtain  asymptotic  expressions  lor  the  electromagnetic  field  reflected  frOTi 
a  surface  of  arbitrary  shape  in  the  case  of  short  waves  (i.e.   for  large  values 
of  the  wave  number  k  =  2n/X),     An  analogous  problem  was  first  considered  by 
V.  A.    Fok  in  19U8.     By  combining  the  Fresnel     laws  of  reflection  and  the  laws 
of  ray  optics  Fok  obtained  formulas  for  tlie  computation  of  the  reflected  field 
of  an  incident  plane  wave  in  the  approximation  of  geometrical  optics  ^  -^ .     A 
generalization  of  these  formtilas  covering  the  case  of  an  arbitrary  incident  wave 
is  found  in  the  paper   [2j  by  V,   A.  Fok,  which  is  quoted  below. 

If  a  diffracting  body  is  an  ideal  conductor,  then  the  values  of  the 
vectors  E(p)   and  H(p)  of  the  field  at  a  point  p  in  space  and  the  current  density 
j  induced  on  the  surface  of  the  body  by 
are  connected  by  the  following  relations: 

ikR 
[  1,  grad  -^ 

S 


(15)  ^(P)  .  ^^(p)  -  ^    /^  rotp  J  [  t,   grad  V  ]  ^' 

S 

(16)  !t(p)  -  !?^(p)  +  ji-     r  [  ^,   grad  ^  ]  dS. 

S 

The  results  of  this  section  are  due  to  U,  K.  Muravev. 


Hi  - 


For  sufficiently  short  wave  lengths  the  current  density  j  is  defined 
up  to  terms  which  decrease  together  with  decreasing  wave  length  (increasing  k) 
by  the  following  relations: 

(17)  i  '   2[n  lt^]     (on  the  lighted  side  of  the  body), 

(18)  J  "  0     (on  the  shaded  side  of  the  body). 

shall  assime  tliat  the  in( 

,   — ^  ikr     .   — ^  ikr      — ^    / —  r  v  — *  1 

(19)  E^-E'V   "o-H'V  •^"'"°=  /1Fo^°]. 

where  E  and  H  characterize  the  amplitude  and  polarization  of  the  incident  wa:v 

on  a  unit  sphere  about  the  source  at  0,  and  r  is  the  xinit  radius  vector  from  0 

to  the  point  in  which  we  define  the  field  ?  ,  H  . 

o'      0 

Using  (I7i   (18),  and  (19)  and  ignoring  terrra  of  order k  «  tte    ve 

(kR)2    ^^ 

obtain  in  place  of  (1$)  and  (16)  the  following  expressions: 

(^,   t,p,  .  k;(p,  .  ^  /I  J  {[^i?  ]  .  l(tSt  ?]))    ^  as. 


^ight 


BQ 


If  the  half -shaded  part  of  the  body  has  small  curvature,  then  for  the  current 
density  in  that  part  of  the  body  we  can  use  the  function  defined  by  V.A,  Fok'-  -' , 
We  restrict  ourselves  to  the  case  when  the  radii  of  curvature  of  the  body  in 
the  half-shaded  part  of  the  body  are  small  in  comparison  with  the  wave  length 
and  we  ignore  the  boundary  effect. 


15- 

ik(r+R) 


(21)  H(p)  .  H^(p)  .  g     J      \tjn  If]]  ^ dS. 

S 

lighted 

•r 
body  to  the  observation  point  p, 

¥oT  an  asymptotic  repi^sentation  of  the  integrals  in  (?o)  and  (2l) 
for  large  values  of  k  we  use  fonmila  (lU). 

Let   the  surface  be  characterized  by  Gaussian  coordinates  u,y  and  lot 
it  be  defined  by  the  equation  v  =»  r(u,v).     The  function  0  defined  above  has, 
in  the  integrals  under  consideration,   the  following  form  (notations  explained  in 
Fig.   2): 


0  -  r  +  R  =  r  +  Yr^  +  r^  -  2(r  r-^), 


Here  r,   denotes  the  distance  between  the  source  and  the  observation  point 
(r,is  a  constant  vector)  and  the  condition  grad  0       [translator:     here  the 
condition  imposed  on  grad  0  is  missing]  yields  the  following  s7/stem  of  equations  for 
the  deterriiiiiation  of  the  coordinates  of  the  stationary  points  on  the  surface  of 
the  body: 

^  St     /-e    ar\ 


^  _  dr  ^  ^  ev  "  V^l  dvj 


which  can  be  written  in  the  form 


16  - 


The  solution  of  this  system  of  equations  yields  the  following  relations 
-I     ^ 
between  the  unit  vectors  R  and  r  (cf.  Fig,  2)  at  the  stationary  point: 


(SU)      Ro  "  "-o  -  2»(°  '■o'' 


(«)      %  ■  ?„, 


where  n  is  normal  to  the  surface. 

Thus,  it  is  possible  to  have  two  types  of  stationary  points: 
Type  1:  the  stationary  point  is  defined  by  equation  (2li)  which  corresponds  to 
the  case  when,  at  an  observation  point  belonging  to  the  lighted  part  of  the 
surface,  the  angle  of  incidence  is  equal  to  the  angle  of  reflection,  and 
Type  2:  the  stationary  point  is  defined  by  condition  (25)  and  is  located  at 
the  intersection  of  the  surface  with  the  line  connecting  the  source  of  the 
field  with  the  observation  point  p  in  the  shaded  part  of  the  space. 

We  assume  that  the  surface  is  such  that  if  p  is  an  observation  point 
belonging  to  the  lighted  part  of  space,  then  there  corresponds  to  it  only  one 
stationary  point  of  the  first  type;  likewise,  if  p  is  an  observation  point 
belonging  to  the  shaded  part  of  space,  then  there  belongs  to  it  only  one 
stationary  point  of  the  second  type. 

Under  these  circumstances  it  is  possible,  using  (lU)j  to  write  (20) 
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and  (21)   as  follows: 


r-     fnH°]    -  R  (t  [nTP])     ^ik(iM-R) 
(26)       E(p)  -  %{v)   -  f\    ^-^ ""-^ ^-=^ H(u,v) 


/ 


V3-? 


0(i), 


St 


(27)   Tt(p)  «=  tip) 


/ 


H(u,t) 


2 


*0(^). 


St 


We  now  express  the  coefficient  /a, a--a_  in  terms  of  the  values  of  the 
local  paraneters  of  the  surface  at  the  stationary  point: 


2^  ,2  I  $)    *     r  ^    -   1  r. 


(28)       a,. 4  I       .4.1s UE 


i^).tti(il]I 


R 


St 


(29)      .,  .    il. 

^        Sudv 


St 


'^r     ^    I  8u     dv/      ^Budv        V'l  Suav/  _ 


a_ 

3uSv 


b3 


St 


(30) 


vS 


St 


.(|)li&ii!J).il4 


^ii)j 


St 


Here  we  assume  that  the  function  9^  -  v  ♦  R  does  not  remain  constant  along  the  contour 
which  divides  the  shaded  and  the  lighted  portions  of  the  surface  of  the  body,  i.e., 
the  sum  of  the  distances  from  the  source  to  the  edge  of  the  body  and  from  the  edge 
*o  the  observation  point  is  not,  for  a  given  form  of  the  body,  a  constant. 
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For  a  stationary  point  of  type  I  the  use  of  condition  (2U)  leads  to  the 
following  expressions  for  the  coefficients  a, ,  a«,  and  a-. 

21 


^2-(l4)h(ii)]-(-o)»'. 


where  E  =  f  |^  j   >  ^  =  f  ~  |^  1  ,  G  «=  f  |^  j   are  the  coefficients  of  the  I 

quadratic  form  of  the  surface  and  D  =  f  n  — |  j,  B     =  (n  .  ^J  ,  D  =  H:  .2-1  1 

\       dm    '  ^    dv    ' 

are  the  coefficients  of  the  II  quadratic  form  of  the  surface. 
Hence 


Here 


+  GD  -  2FD         /I         1    ,     .      ., 

'  'A-  "    "S"  +  o~  1    IS  the  mean  curvature 

g2  lllj_       Rg 


DD  +  D 


'2 


R^ 
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is  the  Gaussian  cunrature. 


It  is  possible  to  show  that  the  expressions 


nii)^° 


(if--( 


h ^ 


are  invariant  with  respect  to  any  transformation  of  the  coordinates  of  the 
surface.  We  transform  the  coordinates  u,t  to  coordinates  ?,>]  in  such  a 
manner  that  at  the  stationary  point  the  curve  >i  =  0  is  the  trace  of  the  inter- 
section  of  the  surface  with  the  plane  of  the  vectors  n  and  r  and  the  curve  ?  =  0 
is  orthogonal  to  the  curve  h  =  0.  Then 


Ai 


!_  «  sin^ 


normal  section  of  the  surface  by  means  of  the  plane  of  incidence. 
V/e  thus  have  for  the  coefficient  in  question: 


/a,. 


M' 


H  <  cos 


^. 


2r  R  cos  e 

r  +  R 


(i^i) 


2o  ^  sin  e 
cos  9  +  -^ 


op   p 

Ur  R  cos  e  1 


1/2 


*  (r  +  R)^   ^2! 
In  the  case  of  a  stationary  point  of  type  II  we  can,  using  condition  (25), 
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reduce  the  expressions  (28)-(30)  to  the  foim: 

^.=(M)h(iii)]. 

whence 


(32)    yL^a^-a^  =/|  +  ^Ucose. 


Thus,  using  relations  (2ii),  (26),  (2?),  and  (3l),  we  can  represent  the 
field  in  the  lighted  part  of  space  as  follows: 


(33)    E(p)  -  E^(p)  -  p  .  2n(n  ?)]^^  B  ^^  *  0  {\ 


(3U)    tiv)   =  H^(P)  +  [?  -  2n(7t  H°)]g^  B  ^  J^       *    0  (|) 


where      r  -.  ,/2 

2  2      1--^/^ 


(3.)   3...i.^[(|^.|^).^j.^3^ 


The  second  terms  in  (33)  and  (3U)  determine  the  aitqiLitude,  polarization, 
and  phase  of  the  reflected  field  in  the  approximation  of  geometrical  optics. 

Using  condition  (25)  and  the  corresponding  formula  (32)  in  computing  (26) 
and  (27)  we  can  show  that 
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^,(P) 


st 


'    /  ?S 


H 


Hjp), 


St 


i.e.,   in  the  approximation  of  geometrical  optics  the  field  in  the  shaded  part 
of  space  is  equal  to  zero. 

In  the  above  considerations  we  assumed  that  the  function  ^  =  r  +  R  is  not 
constant  along  the  boundary  of  the  lighted  and  shaded  domains.     If  r  +  R  -  const, 
on   \      (this  happens,   for  instance,   if  the  observation  point  and  the  source  are 
on  the  axis  of  the  optical  system  which  is  a  mirror  or  a  lens  whose  foim  is 
that  of  a  solid  of  revolution)  then  we  have  the  following  representations  for 
the  field: 

For  the  lighted  part  of  space 

for  the  shaded  part  of  space 
E(P)=S,,„..    0(|), 


where  E  is  defined  by  (33)  and 


This  shows  that  in  this  case  diffraction  phenomena  will   take  place  for  arbitrarily 
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short  wave  lengths  and  E^,  E^^^^,  and  E^^^^  are  all  of  the  sane  order  of 
magnitude . 
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Box  91,  Gardena,   California 
Attn:  Mr.  S.S.  Oleesky 

Library,  Geophysical  Institute  of  the 
University  of  Alaska 
College,   Alaska 

University  of  California 

Berkeley  h,  California 

Attn:  3.  Silver,   Prof.  Engr.  Soienoe 

Division  of  Elec.  Engr. 

Electronics  Research  Lab. 

Electronics  Research  Laboratory 
Unlversitv  of  California 
332  Cory  Hall,    Berkeley  h,   Calif. 
Attn:   J.R.  Whinnery,   Prof. 

California  Institute  of  Technology 
Jet  Propulsion  Laborator;.' 
U800  Oak  Grove  Drive 
Pasadena,   California 
Attn:  I.E.   Newlan 

California  Institute  of  Technology 
1201  E.  California  Street 
Pasadena,   California 

Attn:   Dr.   Charles  H.   Papas 

Carnegie  Institute  of  Technology 
Schenley  Park 

Pittsburgh  13,  Pennsylvania 
Attn:  Albert  E.  Heins 

Cornell  University 

School  of  Electrical  Engineering 

Ithaca,   New  York 

Attn:  Prof.  G.C.  Dalman 

University  of  Florida 
College  of   Engineering 
Gainesville,   Florida 
Attn:   Prof.  M.H.  Latour, 

Engineering  Sciences  Lib. 

Georgia  Institute  of  Technology 
Engineering  Experiment  Station 
Atlanta,   Georgia 
Attn:  Mrs.   J.H.  Crosland 

Georgia  Tech.  Library 

Harvard  University 

Technical  Reports  Collection 

Gordon  McKay  Library,   303A  Pierce  Hall 

Oxford  Street,   Cambridge  38,  Mass. 

Attn:  Mrs.   E.L.   Hufschmidt,   Ub. 

Harvard  College  Observatory 

60  Garden  Street 

Cambridge,  Mass. 

Attnt  Dr.  Fred  L.  Whipple 

UniveriUy  of  IlUnois 
Serials  Dept.  -  2203  Ubrary 
Urbana,   Illinois 

University  of  Illinois 
College  of  Engineering 
Drbana,   Illinois 
Attn:  Dr.   P.E.  Moyes 

Dept.  of  Electrical  Enijr. 


sity 


The  Johns  Hopkins  Univ 
Homewood  Campus 
Department  of  Physics 
Baltimore  18,  Maryland 
Attn:  Professor  Donald  E.   Kerr 

The  Johns  Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,  Maryland 
Attn:  Mr.  George  L.   Seielstad, 
Applied  Physics  Lab. 


Kusichusetta  Institute  of  Technology 
Research  Laboratory  of  Electronics 
Document  Rooai  20B-221 
Canbrldge  39,  Has?. 

Kaasachuaetts  Institute  of  Tech. 

Lincoln  Laboratory 

P.O.  Box  73 

Lexington  73,  Mass. 

Attnt   Henr>-  Straus,   Document  Rb.   A-229 

University  of  Mlchitsn 
Electronic  Defense  Group 
Engineering  Research   Institute 
Ann  Arbor,  Michigan 
Attn:     J. A.  Boyd,   Supervisor 

The  Dniverslty  of  Michigan 
'•lllow  Run  Laboratories 
V.lllo¥  Run  Airport 
Tpsilanti,  Hichigan 
Attni   K.K.  Siegel,   Head 

Theory  and  Analysis  Dept. 

University  of  Michigan 
Willow  Run  Laboratories 
Englneerlnc  Research  Institute 
Willow  Run  Airport 
Tpsilanti,  Michigan 
Attn!  Ubrarlan 

The  Dniverslty  of  Minnesota 

Minneapolis  lU,  Minnesota 

Attn:  Mr.  Robert  H.   Stunra,   Ubrary 

Northwestern  University 
Microwave  Laboratories 
Evanston,   Illinois 
Attn:  Frafassor  R.E.  Bean 

Ohio  State  Dniverslty 
Research  Foundation 
Colunbus  Ifl,   Om.0 
Attn:  Dr.   T.E.  Tlce 

Dept.  of  Electrical  wigr. 

The  University  of  Oklahoma 

Norman,   Oklahoma 

Attn:   Prof.  C.L.  Farrar 

Chairman  Electrical  Engr. 

Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
55  Johnson  Street 
Brooklyn  1,   New  York 
Attn:   Dr.  A. A.  Oliner 

Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
55  Johnson  Street 
Brooklyn  1,   New  Tork 
Attn:     Mr.   A, 


Syracuse  University 
Research  Institute 
Collendale  D-6 
Syracuse  10,   New  Tork 
Attn:  Dr.   C.S.  Grove,   Jr. 

Director  of  Engr.  Research 

The  University  of   I'exas 
Electrical  Engr.  Research  Lab. 
F.%  Box  8026,   University  Station 
Austin  12,   Texas 
Attn:  J.R.  Gerhardt,   Assistant  Dir. 

The  University  of  Texas 

Defense  Research  Lab. 

Austin  12,   Texas 

Attn:  Claude  W.  Horton,   Phys.  Ub. 

The  Dniverslty  of  Toronto 
Department  of  Electrical  Engineering 
Toronto,  Canada 
Attn:  Prof.   G.   Sinclair 

LTI  Eesearch  Foundation 

P.O.  Box  709 

Lowell,  Mass. 

Attn:  Dr.  Charles  R.  Mingins 

University  of  Waahincton 

Dept.  of  Electrical  Engineering 

Seattle  5,  Washington 

Attn:  G.  Held,   Associate  Professor 


Electronicar  ti  vision 
Rand  Corporation 
1700  Main  Street 
Santa  Monica,   California 
Attn:  Dr.  Robert  Kalaba 

National  Bureau  of  Standards 

Waslxinglon,    D.C. 

Attn:   Dr.  W.K.   Saunders 

Applied  Mathonatics  and  Statistics  Lab. 
Stanford  University 
Stanford,   California 
Attn:  Dr.   Albert  H.  Bowker 

Dept.  of  Physics  and  Astronomy 
Michigan  State  University 
East  Lansing,  Hichigan 
Attn:  Dr.   A.  Leitner 

lonsphere  Research  Laboratory 
Pennsylvania  State  College 
State  College,   Pennsylvania 
Attn:  Prof.  A.H.  Waynick,   Director 

Institute  of  Mathanatlcal  Sciences 

25  Waverlv  Place 

New  Tork  3,   New  York 

Attn:  Mrs.   JoAn  Slegal,   Librarian 

School  of  Electrical  Engineering 

Purdue  University 

Lafayette,  Indiana 

Attn:  Professor  K.V.  Schulti 

Univeraitv  of  Tennessee 
Knoxville,   Tennessee 
Attnt  Dr.  Fred  A.   Flcken 

California  Institute  of  Tech. 
1201  E.  California  Street 
Pasadena,   California 
Attn:  Dr.  A.   Erdelyi 

Wayne  University 
Detroit,  Michigan 
Attn:  Professor  A.F.   Stevenson 

Mathematics  Department 
Stanford  University 
Stanford,   California 
Attn:   Dr.  Harold  Levine 

University  of  Minnesota 
Minneapolis  lU,  Minnesota 
Attn:  Prof.  Paul  C.  Rosenbloon 

Department  of  Mathematics 
University  of  Pennsylvania 
Philadelphia  h,    Pennsylvania 
Attn:  Professor  Bernard  Epstein 

Applied  Physics   Laboratory 
The  Johns  Hopkins  University 
6621  Georgia  Avenue 
Silver  Sprin„,  Maryland 
Attn:  Dr.  B.S.  Gourary 

(2)Exchange  and  Gift  Division 
The  Library  of  Congress 
Washington  25,   O.C. 

Electrical  Engineerin  ,  Dept. 
Massachusetts  Institute  of  Technology 
Cambridge  39,  Massachuset-.s 
Attn:  Dr.  L.J.  Chu 

Nuclear  Development  Associates,  Inc. 

5  New  Street 

White  Plains,  New  York 

Attn:  Library 

Lebanon  Valley  College 
Annvllle,  Pennsylvania 
Attn:  Prof.  B.l'..   Bissinger 

Dept.   of  Physics 
University  of  Pittsburgh 
Thaw  Hall 
Pittsburgh  13,  Pa. 
Attn:  Dr.  Edward  Gerjuoy 

Dept.  of  Physics 
Amherst  College 
Amherst,  Massachusetts 
Attn:  Dr.  Arnold  Arons 


California  Institute  of  Tech. 
Electrical  Engineering 
Pasadena,   California 
Attn:  Dr.   Zohrab  A.   Kaprlellan 

Dr.  Rodman  Doll 
209  A  Ermet  Street 
Ypsllantl,  Michigan 

California  Institute  of  Technology 
Pasadena  li,   California 
Attn:  Mr.  Calvin  Wilcox 

Mr.  Robert  Hrockhurst 

Woods  Hole  Oceanographlc  Institute 

Woods  Hole,  Massachusetts 


National  Bureau  uf    • 
Boulder,   Colorado 
Attn:  Dr.  R.   Gallet 


Dr.   Solomon  L.  Schwebel 

36fi9  Louis  Road 

Palo  Alto,    California 

Eastern  Technical  Associates 
38  Main  Street 

Hastlngs-Cn-Hudson,   New  Tork 
Attn:  Dr.  Bernard  Lippnann 

UnlTerslty  of  Minnesota 
The  University  Library 
Minneapolis  lii,  Minnesota 
Attn:  Exchange  Division 

Professor  Bernard  Friedman 

55  HiUtop  Avenue 

New  Rochelle,  New  York 

Lincoln  Laboratory 

Massachusetts  Institute  of  Tech. 

P.O.  Box  73 

Lexington  73,  Massachusetts 

Attn:  Dr.  Shou  Chin  Wang,  Rn.  C-3S1 

Melpar,   Inc. 

3000  Arlington  Boulevard 

Falls  Church,   Vlr,,-inia 

Attn:  Mr.   K.S.  Kelleher,   Section  Head 

Antenna  Laboratory,   EM 

Hq.   Air  Force  Cambridge  Res.  Center 

Laurence  0.  Hanscom  Field 

Bedford,  Massachusetts 

Attn:  Nelson  A.  Logan 

Electronics  Research  Directorate 
Hq.  Air  Force  Cambridge  Res.  Center 
Laurence  G.  Hanscom  Field 
Bedford,  Massachusetts 
Attn:  Dr.   Philip  Newman,   CRRK 


Crosley  ATCO  Res,  Labs. 
257  Crescent  Street 
Waltham,  Massachusetts 
Attn:  Mr.  "l.C.  Serson 

Hoffman  Labs.,  Inc. 
Advanced  Development  Section 
3761  South  HiU  Street 
Los  Angeles  7,   Calif. 
Attn:  Dr.  Richard  B.   Barrar 

Columbia  University 
Hudson  Laboratories 
P.O.  Box  239 
11j5  Palisade  Street 
Dobbs  Ferry,   New  Tork 
Attn:    Dr.  H.W.   Johnson 

Institute  of  Fluid  Dynamics 

and  Applied  Math,    Univ.  of  Maryland 

College  Park,  Maryland 

Attn:  Dr.   Elliott  Hontroll 


Dept.  of  Electrical  Engineering 
Washington  University 
Saint  Loviis  5,  Missouri 
Attn:  Prof.   J.   Van  nladel 

Dept.   of  the  Navy- 
Office  of  Naval  Res.  Branch  Of  fie 
1030  E.  Green  Street 
Pasadena  1,   California 

Brandeis  University 
Waltham,  Massachusetts 
Attn:  Library 

General  Electric  Company 
Microwave  Laboratory 
Electronics  Division 
Stanford  Industrial  Park 
Palo  Alto,   California 
Attn:  L-lhrary 


Smyth  Research  Associates 
3930     Uth  Avenue 
San  Diego  3,  California 
Attn:  Dr.   John  B.   Smyth 

Electrical  Engineering 
California  Inst,   of  Technology 
Pasadena,    California 
Attn;  Dr.  Georges  C.  Weill 

Naval  Research  Laboratory 
Washington  25,  D.C. 
Attn:  Dr.  Henry  J.  Passerini 
Code  5278  A 

Dr.  George  Kear 

10585  N.   Stelllng  Road 

Cupertino,  California 

Brooklyn  Polytechnic 

85  Livingston  Street 

Brooklyn,   Hew  York 

Attn;  Prof.  Nathan  Marcuvitz 

Dept.   of  Electrical  Engineering 
Brooklyn  Polytechnic 
85  Livingston  Street 
Brooklyn,   New  York 
Attn:  Dr.  Jerry  Shmoys 

Dept.   of  Mathematics 
University  of  New  Mexico 
Albuquerque,   New  Mexico 
Attn:  Dr.  I.   Kolodner 

W.L.  Maxon 

1^60  West  3Uth  Street 

Hew  York,   N.Y. 

Attn;  Dr.  Harry  Hochstadt 


